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Abstract 



The physics of the fractional quantum Hah effect is the physics of interact- 
ing electrons confined to a macroscopically degenerate Landau level. In this 
Chapter we discuss the theory of the quantum Hall effect in systems where 
the electrons have degrees of freedom in addition to the two-dimensional or- 
bital degree of freedom. We will be primarily interested in the situation where 
a finite number of states, most-commonly two, are available for each orbital 
state within a degenerate Landau level and will refer to these systems as multi- 
component systems. Physical realizations of the additional degree of freedom 
include the electron spin, the valley index in multi-valley semiconductors, 
and the layer index in multiple-quantum-well systems. The consideration of 
multi-component systems expands the taxonomy of incompressible states and 
fractionally charged excitations, and for example, leads to the appearance of 
fractions with even denominators. More interestingly, it also leads us to new 
physics, including novel spontaneously broken symmetries and in some cases, 
finite temperature phase transitions. We present an introduction to this rich 
subject. 



To appear in: Novel Quantum Liquids in Low- Dimensional Semiconductor Structures^ edited 
by Sankar Das Sarma and Aron Pinczuk (Wiley, New York, 1995). 
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I. INTRODUCTION 



The fractional quantum Hall effect ^""^ is a remarkable example of strong correlations in a 
two-dimensional electron gas (2DEG). In zero magnetic field, a dimensionless measure of the 
strength of the Coulomb interaction for a system with dielectric constant e, Fermi energy ep 
and Fermi wavevector kp is 

A = ^^, (1) 

which is small in the limit of high density. In this limit one can frequently treat the effects 
of the Coulomb interaction perturbatively. Physically this can be visualized as being due to 
the electrons moving rapidly past each other at the Fermi velocity and thus not scattering 
so strongly as they would at lower densities. A strong magnetic field completely changes 
this situation. Semiclassically, the rapid motion of the electrons is converted into circular 
cyclotron orbits. The particles now scatter strongly from each other and in fact semiclas- 
sically do not move except under the E x B drift induced by their mutual interactions. A 
full quantum treatment of the motion shows that the kinetic energy is quenched and now 
occurs only in discrete values (n -|- l/2)hu!c, where n is the Landau level index. Since the 
kinetic energy within a given Landau level is completely degenerate, the Coulomb interac- 
tion inevitably induces highly non-perturbative effects. (The Landau level degeneracy is 
— BA/^o where $o = hc/e is the magnetic flux quantum, and A is the area of the 
system.) 

The essential feature of the fractional quantum Hall effect is a condensation of the elec- 
trons into special highly correlated states^ which minimize the Coulomb energy by having 
the electrons avoid each other as much as possible. These states are characterized by an 
unusual topological order^^ which costs a flnite amount of energy to break. Hence the fluid 
is incompressible and has an excitation gap both for its charged^ and neutraP*^' excitations. 
P. W. Anderson has characterized this state as a Mott insulator induced by the magnetic 
field. 12 

The essence of this phenomenon is captured in a remarkable class of wave functions first 
constructed by R. B. Laughlin 

N f 1 ^ 1 

i^rr^iZi, Z^, Z^) = l[iZ^- Z,r exp -- E \Zk? . (2) 

Here, since we are in two dimensions, we are using the dimensionless complex number Z — 
(x+iy) /ito represent the position vector {x, y) in units of the magnetic length I — ^Jhc/eB. 
This wave function describes spinless fermions in the lowest Landau level (in the symmetric 
gauge). To satisfy the analyticity requirement placed on the wave function by the constraint 
of being in the lowest Landau level, the parameter m must be an integer. To satisfy 
the antisymmetry requirement for fermions, m must be odd. Laughlin's plasma analogy^ 
shows that the parameter m fixes the Landau level fiUing factor to he u = N/N^ — 1/m. 
Experiments^ has indeed observed gapped quantum Hall states at filling factors ly — 1, 1/3, 
and 1/5. 
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It is clear that Laughlin's wave function builds in good correlations because it vanishes 
as |rj — Tjl"^ when any two particles i and j approach each other. Thus there is only a 
small amplitude for the particles to be near each other, and the Coulomb energy is lowered. 
Note that no pair of particles ever has relative angular momentum less than m. Hence the 
Laughlin function is a zero-energy exact eigenfunction for Hamiltonians with the appropriate 
finite number of non-zero Haldane pseudopotentials,^^ Vn'- 

H^J2ynT.Pn[t,j], (3) 
n<m i<,j 

where Pn[iij] is the lowest Landau level (LLL) projection operator for the relative angular 
momentum state n of particles i and j. In the lowest Landau level, relative angular momen- 
tum is proportional to the square of the separation. Hence the Laughlin wave function is 
very nearly an exact ground state for any sufficiently short-range repulsive interaction. 

In addition to the primary filling fractions v = 1/m, numerous other fractions have 
been observed, all of which (for single component systems) have odd denominators (again 
because of the Pauli principle) . These phases have been explained in terms of a hierarchical 
picture using bosonic,^^ anyonic,^^ and fermionic^' representations. More recently Jain^^'^^ 
has discussed an appealing composite fermion picture. N. Read has argued that all of 
these representations are mathematically equivalent^'' and contain the same physics. Which 
representation is most convenient depends on circumstances. Jain's approach has inspired 
several new and important experiments.^^ 

Our purpose here is to consider the nature of the various phases which can occur in multi- 
component systems. There are several physically different realizations of systems with extra 
degrees of freedom which require a multi-component representation, and important early 
work on this problem was done by Halperin.^^'^^ 

The first and simplest example is that of ordinary electron spin. In free space for electrons 
with g factor 2 (i.e., neglecting QED corrections) the Zeeman sphtting q/jLeB is precisely 
equal to the Landau level splitting huc. If this were true in quantum Hall samples, then 
even a non-interacting system with filling factor u = 1 would have a large excitation gap 
for fiipping spins and the ground state would be fully polarized at low temperatures. In 
this case spin excitations are frozen out and we can treat the electrons as being effectively 
spinless. However in the solid state environment of the 2DEG, two factors conspire to to 
make the effective g-factor much smaller in many semiconductors and particularly in GaAs 
samples in which almost all fractional quantum Hall studies have been done. The first is that 
the small effective mass (m* ~ 0.068 in GaAs) increases the cyclotron energy by a factor 
of approximately 15. Secondly, spin-orbit coupling causes the spins to tumble and reduces 
their coupling to the external magnetic field by roughly a factor of 4. Thus the ratio of the 
Zeeman splitting to cyclotron splitting is reduced from unity to about 0.02 in GaAs. The 
spin-orbit contribution is pressure dependent and may also be further reduced as a result of 
size quantization effects in narrow quantum wells. 

For small enough g factor and weak enough magnetic fields, spin fluctuations become an 
important dynamical degree of freedom and we must use a two-component wave function to 
describe the system. While the Coulomb force is spin-independent, we shall see below that 
exchange effects can lead to spontaneous ferromagnetism as well as to gapped 'local singlet' 
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spin liquids, which are, loosely speaking, itinerant antiferromagnets. [Unlike antiferromag- 
nets however, they do not break translation symmetry nor do they have gapless Goldstone 
modes). Ferromagnetism turns out to be important at filling factor v = 1 where Coulomb 
exchange effects are much more important than originally realized. 

A different class of quantum Hall states with gaps has been observed experimentally^^"^^ 
at filling fractions with even denominators such as = 5/2. These have been argued^' ^^"^^ to 
be either manifestations of spin liquid states or special 'p-wave pairing' spin-polarized states 
known as pfaffian states. One useful signature of spin effects in 2DEG's is that they are 
sensitive to tilts of the magnetic field. To a first approximation, the orbital degrees of freedom 
are sensitive only to the perpendicular component of the field because of the 2D confinement, 
while the Zeeman splitting is proportional to the total magnetic field. (However, the coupling 
of orbital degrees of freedom to the parallel component of the magnetic field is not completely 
negligible in typical quantum wells^^ and this can considerably complicate the situation.) 

A second example of a multi-component system is found in silicon where the conduction 
band minimum occurs not at the F point (the zone center) but rather at six symmetry 
equivalent points lying near the zone boundary along the principal cubic directions. Thus 
electrons doped into the conduction band of Si must be described by a six-component wave 
function (if we ignore spin). The presence of the oxide barrier in a Si MOSFET device and 
the enormous electric field perpendicular to it (which is used to confine the electrons into 
the inversion layer) breaks the cubic symmetry and lowers the energy of two of the valleys. 
For typical electron densities, only these two valleys are occupied, thus yielding a system 
which is effectively two-component and has SU(2) symmetry just like a spin-1/2 system.^'' 

A third example which will be extensively discussed here and in the Chapter in this 
volume by J. P. Eisenstein^^ occurs in double quantum well structures. With modern MBE 
techniques it is possible to grow GaAs heterostructures containing two 2DEG's separated by 
a distance comparable to the spacing between electrons within each layer. Remarkably, it is 
also possible to make separate electrical contact to each layer. A closely related system is a 
single wide quantum well in which the two lowest electric subbands are nearly degenerate.^" 
We will make the (only approximately correct) assumption throughout our discussion that 
the low energy physics of a single wide well can be mapped onto that of a double well 
with appropriately chosen parameters. In all two component systems it is useful to define a 
pseudospin representation in which spin up and down refer to the two possible values of the 
layer index (or subband index) for each electron. We will frequently frame our discussion 
in a spin or pseudospin language; the reader should be aware that the discussion applies 
equally well to double-layer systems. When the distinction between the U(l) and SU(2) 
symmetries of the interaction term in the Hamiltonian is important, we will say so. 

Because the intra- and inter-layer Coulomb matrix elements are different, double well 
systems do not have full SU(2) symmetry, but rather only U(l) symmetry associated with the 
conservation of the charge difference between the two layers (assuming there is no interlayer 
tunneling). ^^"^^ These systems exhibit gapped quantum Hall states with both even and 
odd denominators. In addition, there are gapless XY ordered phases. These phases are 
destroyed above some critical temperature by a Kosterlitz-Thouless phase transition. This 
phase transition (not yet observed experimentally) is the first example of a finite temperature 
phase transition in a quantum Hall system. 

This Chapter is organized as follows. In Section || we introduce the physics of incom- 
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pressible states in multi-component cases by presenting the appropriate generalizations of 
Lauglilin's many-particle wavef unctions. Section |TTl| briefly directs the reader to references 
on the Chern Simons effective field theory approach to these problems. Section |V| discusses 
fractional charges in multicomponent systems. Section ^ discusses collective modes from 
the point of view of the single-mode approximation. In both these sections we will see that 
some of the multi-component wavefunctions have correlation functions whose qualitative 
properties deviate from the norm. These differences lead to changes in the nature of the 
fractionally charged excitations and in the collective modes. The qualitative difference is 
associated with a broken symmetry which occurs in some cases as we discuss in Section |V^. 
In our view, it is in the properties of these broken symmetry states that some of the most 
interesting new physics in mult i- component quantum Hall states is revealed. Section |VI1 
discusses the field-theoretic gradient expansion approach to the broken symmetry case. In 
double-layer systems, the broken symmetry state has spontaneous phase coherence between 
the electrons in different layers even when these layers are isolated apart from inter-layer 
Coulomb interactions. Section |^ deals with the external symmetry breaking introduced 
by interlayer tunneling. Spontaneous coherence in double-layer systems leads to remarkable 
effects upon tilting the magnetic field away from the normal in a double layer system. Some 
of these effects are discussed in Section ^. Finally Section ^ presents a summary of the 
central ideas. 



II. MULTI-COMPONENT WAVE FUNCTIONS 

In this section we discuss wavefunctions for spin or pseudospin 1/2 particles which can 
be written in the symmetric gauge in the form 

^[Z : x] = AmZ]a, ■ ■ ■ a^v^An " " " A^j]' (4) 

where A is the antisymmetrization operator, [i] = N-^ + i, and and jSk are the spinors 
for the kth electron aligned respectively parallel and antiparallel to the Zeeman field (which 
we take to be in the z direction). Thus ^[Z] is the orbital wave function for the spin 
configuration in which the first N-^ electrons have spin "up" and the remaining electrons 
have spin "down".^^ 

The two-component orbital wave functions originally proposed by Halperin^^'^^ have a 
form analogous to the Laughlin functions 

^m,m',n[Z] = H ~ ^i)"" II ("^W ~ ^[l])"^ 

i<j<N^f k<l<Ny 

xl[l[iZ^-Z[,^rl[exp{-\Zf/A}, (5) 

a=l 6=1 s=l 

where Z^ = (xk + iyk)/^ is the 2D layer coordinate of the kth electron expressed as a complex 
number, and m and m' are odd integers. ^rn,m',n[Z] excludes relative angular momenta less 
than m between up-spins, less than m' between down-spins, and less than n between an 
up-spin and a down-spin. Hence the same arguments presented above for the Laughlin 
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wave function can be used to motivate the idea that the Halperin wave functions are good 
approximations to the ground state for short-range repulsive potentials. 

In the single component case, Laughlin wavefunctions are expected to accurately ap- 
proximate the ground state for m = 3 and m = 5(z/ = l/3 and u = 1/5); at smaller filling 
factors the Wigner crystal state takes over. In the two-component case, the realization of 
a Laughlin state at a particular filling factor may be dependent on other parameters of 
the system such as the layer separation in double-layer systems. One slight complication 
is presented by the SU(2) symmetric case.^^'"^^ Here the total spin commutes with the 
Hamiltonian and one may, without loss of generality, require that the energy eigenstates 
simultaneously be eigenstates of S^. This requirement is not satisfied in general by the 
Halperin wave functions. It is easy to pick out two special cases (among others) that do 
work however. The choice m = m' = n yields a fully antisymmetric spatial wave function 
and hence implies a fully symmetric spin function. This immediately tells us that we have 
a fully aligned ferromagnetic state with total spin quantum number S = N/2. On the 
other hand, the choice {m, m',n} = {1,1,0} corresponds to a simple Slater determinant 
with both spin states of the lowest Landau level fully occupied, giving z/ = 2. Hence it is 
automatically a spin singlet. We can generalize this to {m,m',n} = {m, m, (m — 1)} since 
this corresponds to simply multiplying the filled Landau level function by a fully symmetric 
spatial polynomial. Thus these states are also spin singlets. 

Using the fact that every extended single-particle orbital in the lowest Landau level 
involves a polynomial with a fixed (average) density of zeros^^'^^ given by B/^q (where 
$0 = hc/e is the flux quantum) we may derive a pair of equations for the density p of each 
component 

B 

— = mpi + npi, 

= m'p^ + npp (6) 

^0 

From this we obtain for the filling factors 

m' — n 

mm' — n'^ 
m-n 

n = 2- (7) 

mm' — 

Partial and total filling factors for some of these two-component Laughlin states are listed 
in Table |. At these partial filling factors <^rn,m' .n[Z] is unique in the sense that it is the only 
wave function which excludes its corresponding low relative angular momentum channels 
and, just as in the case of Laughlin states in one-component systems, we may expect that 
these wave functions will be nearly exact groundstates for any sufficiently short-ranged 
repulsive interaction. In the SU(2) invariant case it seems that we should^° however require 
that the wave functions be eigenstates of the total spin operator St and that their Zeeman 
energy be not too unfavorable. There is numerical evidence for instance, that the {3, 3, 2} 
state which has z/ = z/| + 1/^ = 2/5 has a lower energy than the usual hierarchical state if (and 
for typical field strengths, only if) one ignores the Zeeman energy. For further references and 
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detailed discussion on this point and other topics related to spin in the FQHE, the reader 
is directed to the book by Chakraborty and Pietilainen.^ 

We note that Eq.(|^ is ill-defined if mm' — n"^ vanishes as it does, for example, in the fully 
ferromagnetic {m, m, m} states. In this case, however, one can compute the filling factor by 
simply noting that the fully spin aligned system has an orbital wave function equivalent to 
the Laughlin function at total filling factor u = 1/m. The relative filling factors of the two 
components is necessarily ill-defined because of the SU(2) rotational symmetry. There are 
25* + 1 = iV -|- 1 orthogonal but macroscopically degenerate states differing only by their 
quantum number. Thus we have 

1 

n + ^'i = — ' 
m 

- n = J^- (8) 

The degeneracy of these states leads to a broken symmetry which we discuss in greater detail 



in Section VI 



Generalizations of the Halperin states can be made to the case of an arbitrary number 
of components. ^^"^^ This would have application for example to a superlattice of closely 
spaced quantum wells, should these become technologically feasible to produce at some 
point in the future. It is known experimentally^^"^^ that there exists an incompressible (but 
unusually delicate) Hall state at filling factor z/ = 5/2. This state has also been observed 
numerically for various artificially chosen interaction models, but not however, for a pure 
Coulomb interaction.^' ^^'^^'^^ It has been argued that the ground state at this filling factor 
is not fully spin polarized because the Hall state is easily destroyed by tilting the magnetic 
field at constant filling fraction. In principle, the same state should be observed at filling 
z/ = 1/2 because z/ = 5/2 = 2 + l/2 has both spin states of the LLL filled and these electrons 
are essentially inert, leaving an effective filling factor u^g = 1/2 in the next Landau level. 
However the analogous state is not observed at z/ = 1/2 because, to reach this lower filling 
factor (in the same GaAs sample), it is necessary to increase the perpendicular component of 
the magnetic field by a factor of 5. This increases the Zeeman energy advantage of the fully 
spin polarized state (which is gapless and does not exhibit a Hall plateau) and makes it the 
ground state. Thus non-observation of a quantum Hall state at u = 1/2 provides additional 
evidence that the z/ = 5/2 state is not fully spin polarized. As described in the Chapter 
in this volume by Eisenstein, tilted field experiments can also help identify spin-reversed 
quasi-particle excitations above polarized ground states.^' 

Haldane and Rezayi^^ have proposed a two-component spin-singlet wave function to 
explain the existence of the u = 5/2 quantum Hall plateau. This wave function is an exact 
ground state for the so-called 'hollow-core' model^^ and may be written in two different ways 
by modifying two different Halperin states, each of which has filling factor one half.^^ The 
first uses the Halperin fermionic {3, 3, 1} function 

$HC = $331 per|M|, (9) 
where M is an A^/2 x A^/2 matrix whose ij element is given by 

M,, = (Z,-Z[,])"\ (10) 
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The permanent of the matrix, per|M|, is by definition, just hke the determinant except 
that there are no minus signs for odd permutations. The subtle effect of this permanent on 
the wave function is to cause it to be an eigenfunction of total spin (with S = 0) without 
changing the density (at least in the thermodynamic limit). 

Note that the permanent causes some up and down spin particles to have a finite prob- 
ability of having relative angular momentum zero. It turns out that no particles ever have 
relative angular momentum 1 in this state. Hence this wave function is an exact zero energy 
singlet ground state for the 'hollow core' potential model 

Kn = V6m,l. (11) 

Despite the unphysical appearance of this model, it has been argued that it might capture 
the correct physics when the form of the effective pseudopotentials in the second Landau 
level is taken into account. 

A second way to write the same state uses the (222) bosonic Halperin wave function 

$HC = $222Det|M|, (12) 

where M is an N/2 x N/2 matrix whose ij element is given by 

M,, = {Z, - Zyi)'". (13) 

A curious mathematical identity^^ allows one to show that these two representations are 
precisely equivalent. 

It is possible that orbital effects^^ confuse the tilted field test for spin unpolarized states 
and the u = 5/2 is actually spin-polarized. A competing spin-polarized 'pfaffian' state 
developed by Read^ and also studied by Greiter et al.^^ is a kind of 'p-wave' paired state. 
It is not known for certain at this point what the true nature of the very delicate 5/2 state 
jg 38,54 ixiay be one of the proposed states^' ^^"^^ or something completely unknown. 

Jason Ho has recently considered interesting connections between the internal order in 
these types of wave functions and analogous order in superfiuid ^He.^'' In particular he has 
discussed the "incompressible" deformation of the various states into each other, connecting 
for example the {3, 3, 1} state and the pfaffian state.^^'^'''^^ 



III. CHERN SIMONS EFFECTIVE FIELD THEORY 

One interesting approach to the quantum Hall effect in general, and multi-component 
systems in particular, is the Chern Simons effective field theory. Unfortunately space limi- 
tations prevent us from discussing this approach in any detail. The reader is directed to the 
Chapter in this volume by B. I. Halperin, the references therein and to the many refer- 
ences in what is now a vast literature.^^'^^"^^'^^'^""^^ A succinct and introductory summary 
of the bosonic representation for double layer systems is given in Ref . [46] . We present a brief 
discussion of the collective mode predictions of the Chern Simons approach for double-layer 
systems in Section |V| of this Chapter. 
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IV. FRACTIONAL CHARGES IN DOUBLE-LAYER SYSTEMS 



The fractional quantum Hall effect occurs because, at particular filling factors, electrons 
in a partially filled Landau level are able to organize themselves into such strongly correlated 
states that the energy cost of making decoupled particles and holes remains finite, even in 
the thermodynamic limit, i.e., there is a charge gap. It is unusual to have a charge gap 
due entirely to electron-electron interactions (i.e. in a tranlationally invariant continuum 
system) although the example of superconductivity is familiar. In the fractional quantum 
Hall effect, not only do interactions produce a charge gap, but the free-charges responsible for 
the thermally activated dissipation measured experimentally contain only a fraction of the 
charge of an electron. The fact that sharply-defined fractional charges occur in the fractional 
quantum Hall effect can be understood as a necessary consequence of the quantization of the 
Hall conductance.^' (See the related discussion in Section |V11| .) It can also be understood 
in terms of the variational wavefunctions introduced by Laughlin in his pioneering work on 
the theory of the fractional quantum Hall effect.^ In this section we discuss the fractionally 
charged excitations of the '^rn,m',n states by generalizing the plasma arguments made by 
Laughlin for single-layer systems. 

The \l/m,m',n wavefuuctious have the property that pairs of electrons are excluded from 
certain relative angular momentum states. A low-energy charged state must have a localized 
excess or deficiency of charge without destroying the energetically favorable correlations 
associated with the relative angular momentum state exclusions. For a single-component 
system, Laughlin suggested that an accurate approximation to the many-body wavefunction 
for a state with a charged hole at the origin could be obtained simply by multiplying his 
wavefunctions for incompressible states by the factor Hi Zi. For two-component systems this 
argument has an obvious generalization. We can produce two, in general different, charged 
excitations centered on the origin, by multiplying the orbital many-particle wavefunction by 
the product of Zi for all electrons in one pseudospin state. This operation is the variational 
wavefunction equivalent of introducing an unattached flux-tube in Chern-Simons theories 
and the conclusions we reach below can equally well be obtained by using an algebraically 
equivalent argument in that language. The plasma analogy results from writing the quantum 
distribution function, the square of the many-body wavefunction, as a classical statistical 
mechanics distribution function for interacting particles in an external potential 



The classical systems that result are generalized two-dimensional Coulomb plasmas and 
it is convenient in discussing them to adopt the convention of using Roman indices for one 
component of the plasma and Greek indices for the other component. With this notation, 
the trial wavefunctions we consider for charged excitations are 




(14) 



^ m,m' ; 



(15) 



i=l 



and 




(16) 
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These trial wavefunctions clearly reduce the density near the origin without ruining the good 
correlations in '^rn,m',n- At first sight it might appear that only the Roman particle density 
is reduced in ^, „ and only the Greek particle density is reduced in ^, „, but this is 
not the case in general because of the correlation factors. 
The classical potential energy corresponding to \E'^ ^, „ is 



U^^^,^„ = mJ2{-2ln\Z,-Z,\) + m'J2i-2ln\Z^-Zp\)+nY.{-2\n\Z,-Z^\) 

i<j a<f3 i,Q 

This potential is a generalized two-dimensional Coulomb plasma'^^ in the sense that the 
coupling constants outside the sums in the interaction terms are not constrained to be 
the products of charges for the two species, i.e., we allow n ^ (mm')^/^. This difference 
results in long-range interactions in the plasma which depend on the density of each species 
separately, rather than just on the total 'charge'. U^^, ^ is the potential energy function for 
a system consisting of Roman and Greek particles. All particles have repulsive mutual two- 
dimensional Coulomb interactions with coupling constant m between two Roman particles, 
coupling constant m' between two Greek particles, and coupling constant n between a Roman 
particle and a Greek particle. All particles are attracted to a neutralizing background which 
can be considered to have resulted from interaction with unit coupling constant with non- 
responding particles of uniform charge density {2it£'^)~^. For U^, only Roman particles 
interact with unit coupling constant with an impurity particle located at the origin. 

The charge densities induced in each species of particles by the impurity can be calculated 
using the perfect screening properties which result from the long-range interactions of the 
plasma. Far enough from the impurity the direct long-range interaction must vanish for 
each species of particle; i.e., the sum of the impurity charge times its coupling constant plus 
the induced charges in each plasma component times the coupling strength for that plasma 
component must vanish so that 

(17) 

In Eq.(p]7D is the contribution, in units of the magnitude of the electron charge, to the 
quasiparticle charge from Roman particles and Cq is the contribution from Greek particles. 
Eq.(|T7D can be solved for and Cq and the total quasiparticle charge = + Cq: 

mm' — mm' — mm' — 

The fractional charges for \E'^ ^, „ differ only through the interchange of m and m': 




-n n m ^ m — n 



= 4 = — ^- (19) 



mm' — mm — n'^ mm' — 

Fractional charges calculated from these expressions are listed in Table 
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We note that the total charge of what is presumably the lowest energy charged excitation 
at each filling factor, has a value e/q where q is the denominator of the fractional total filling 
factor. When the two- components of the incompressible Hall fiuid are correlated, a reduction 
of charge density in one-component leads to an increase of charge density in the other 
component. The total charge thus tends to consist of partially canceling contributions from 
the two-layers. This cancellation reaches its extreme limit for the case where m = m' = n, for 
which the total charge of the excitation is well defined but its separation into contributions 
from separate components cannot be fixed by the perfect screening requirement on the 
plasma. We will see later that what is behind this behavior is the existence of long-range- 
order in the '^m,m,m wavefunction. This long range order has delivered a bonanza of new 
physics in two-component systems, which will be the focus of much of this Chapter. 



V. COLLECTIVE MODES IN DOUBLE-LAYER QUANTUM HALL SYSTEMS 

In this section we will discuss both intra-Landau-level and inter-Landau-level (cyclotron) 
collective (neutral) excitations of the incompressible ground states whose origin we have 
explained in previous sections. Our discussion is based on the projected single-mode- 
approximation, which has proved extremely useful^^'^^' ''^ in understanding the nature of 
the collective mode structure in single-layer systems. The projected single-mode approxi- 
mation can be appropriate in the strong magnetic field limit where there is little Landau level 
mixing in either the ground state or the low-lying excited states. Many-body eigenstates of 
the system can then be distinguished by the quantized integer number of units of by 
which the kinetic energy exceeds the minimum value NhuJc/2 (where N is the number of 
particles). The single- mode-approximation for the collective energy spectrum follows from 
the assumption that there is a unique many-body state, 1^1/^), with energy EJ^ within each 
quantized kinetic energy manifold, which is coupled to the ground state by the one-body 
density operator: 

Pk = '^exp{ik-ri). (20) 

i 

[We will employ complex number notation [k = + iky) for two-dimensional vectors when 
convenient.] 

To use the single-mode approximation it is necessary to separately evaluate contributions 
to moments of the dynamic structure factor from transitions involving different numbers of 
quantized kinetic energy units. We write for the dynamic structure factor 

s{k,e) = J2sn{k,e), (21) 

n 

where 

Sn{k,e) = ^5^|(^,,„|pfc|^o)r^(e-^^,n + ^o), (22) 

and A is the area of the system. Here |^'i,ri) is an exact eigenstate of the Hamiltonian 
involving n excess quantized kinetic energy units. In the single-mode approximation it is 
assumed that only a single eigenstate contributes to the sum in Eq. (p^) , 
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Sn{k, e) = l|(iI/fc,„|p,|il/o)|25(e - E« + Eo). (23) 

Both the matrix element and the energy which appear in this expression have physical 
significance. The matrix element determines how strongly one-body external probes (for 
example, far infrared or microwave radiation) couple to the collective excitation. The en- 
ergies of the collective modes can be measured in transmission experiments or in inelastic 
light-scattering experiments.^^ The thermodynamics and linear response functions of the 
system depend qualitatively on the presence or absence of collective modes whose energies 
vanish in the limit of long wavelengths. 

Two moments of Sn{k,e) are relatively easy to evaluate and we will use these two mo- 
ments to determine both the matrix element and the collective mode energy: 



Jo 



desn{k,e), (24) 
and 

fn{k)= r deesn{,k,e). (25) 

JO 

Given s„(/c) and fnik), we have for the matrix element 

\{^l\p,\^o)\'^Asn{k), (26) 

and for the excitation energy 

A]^ = E--Eo = Uk)/sr,{k). (27) 

Since ph=Q is a constant, and the ground and excited states must be orthogonal, it follows 
that at long- wavelength s„(A;) must vanish at least as fast as k'^. Long wavelength probes, 
like far infrared radiation, produce^^ observable coupling to a long-wavelength collective 
mode only if Sn{k) cx k'^ at long- wavelengths. We will refer to modes for which Sn{k) ex A;^ as 
dipole active and those for which Sn{k) vanishes with a higher power of /c as dipole inactive. 
The quantity fnik) is proportional to the product of the square of the matrix- element and 
the excitation energy and in analogy with atomic physics we refer to this quantity as the 
projected oscillator strength. We will refer to Sn{k) as the projected static structure factor. 
The usual 'f-sum rule', valid with or without a magnetic field, states that for parabolic bands 
with effective mass m* 

Y,fn{k) = ^h^k^/2m\ (28) 

To evaluate projected 'f-sum rules' we exploit the property that the single-particle Hilbert 
space of a charged particle in a magnetic field can^^ be considered as the product space of a 
factor space in which states are distinguished by the number of quantized kinetic energy units 
in the cyclotron orbit and a factor space for the cyclotron-orbit-center degree of freedom 
which exists within each Landau level and is responsible for the macroscopic Landau level 
degeneracy. To separate the dynamic structure factor into contributions associated with 
different quantized kinetic energies we write the kinetic energy operator in the form 
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P. = Ep^^ (29) 

n' ,n 

where 

pt^'' = Y.\^')^{n\,G''''-{k)B,{k). (30) 

i 

The G" '"(A;) are related to Laguerre polynomials, ''^ the sum over i is over particle labels, and 
Bi(k) is a factor coming from the projection of pk onto a single Landau level and operates 
on the intra-Landau-level degree of freedom of particle i. The commutators which appear 
below are evaluated by using Eq. (^) and the identity, 

Bi{h)Bi{k2) = expiklk2/2)Bi{k, + k^). (31) 

At strong magnetic fields, only p^'° contributes to s„(A;) and fn{k). We will restrict our 
attention to the cases of intra-Landau-level (n = 0) and magnetoplasmon [n = 1) collective 
excitations. (Excitation modes out of the LLL with n > 1 are never dipole active.) 

Using completeness relations of the many-body eigenstates it follows from the above 
definitions that 

so{k) = ^(^o|p°J^ri*o) (32) 

and that 

fo{k) = ^{mp'-lnpn^o) (33) 

where V is the electron-electron interaction term in the Hamiltonian. From general prop- 
erties of (translation invariant) many-body eigenstates within the lowest Landau levels it 
can be shown^° that for small k, so{k) oc k'^ so that the intra-Landau level collective mode 
is dipole inactive. After a somewhat laborious calculation, involving repeated application 
of Eq.(^), the n = projected oscillator strength can be expressed in terms of so{k) and 
one finds the result that fo{k) oc /c*^, independent of any details of the electron-electron 
interaction or the ground state wavefunction, so that the magnetoroton intra-Landau-level 
collective modes have a gap at long-wavelengths: 

limA°^0. (34) 

k—i-O 

On the other hand an elementary calculation based on Eq.(pTD and the strong- field-limit 
assumption that the ground state lies entirely within the lowest-Landau-level subspace of 
the full Hilbert space, implies that 

N k^f^ 

.i(fc) = -^exp(-|H|V2). (35) 

Similarly^^ 

Mk) = huMk) + ^{^o\[[p':lnpl'']\^o), (36) 
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and the second term on the right-hand-side of Eq. (|36D can be shown to vanish as for small 
k. The magnetoplasmon mode is dipole active and, in the long-wavelength limit, completely 
exhausts the full f-sum rule. For a single-layer, the magnetoplasmon mode is the only dipole- 
active mode and its energy is not shifted from by electron-electron interactions. These 
behaviors result from the conservation of particle number and invariance under translation: 
a long wavelength electromagnetic field couples only to the cyclotron motion of the center 
of mass of the system.'''^ 

We are now prepared to discuss how these results are altered in double layer systems. 
We restrict our attention to the case where the two 2DEG's are identical and tunneling 
between them may be neglected. In this case the number of electrons in each layer is a 
good quantum number and collective modes corresponding to the sum and difference of the 
density oscillations in the two layers decouple. To generalize the projected single-mode- 
approximation to the double-layer case we evaluate separately projected oscillator strengths 
for both sum (in phase) and difference (out of phase) modes. We will find that the behavior 
of the sum modes for double-layer systems is similar to the behavior of the modes of a 
single-layer system while the behavior of the difference modes departs from this pattern. 
For the difference modes, both n = and n = 1 modes are dipole active. The n = 1 mode 
is shifted from huc and the n = mode usually has a finite energy as ^ 0. An exception 
occurs for those ground state wavefunctions which have a type of long range order which 
we have not yet made explicit. This long range order is associated with a broken symmetry 
ground state. These broken symmetry ground states are, arguably, responsible for the most 
surprising and appealing new physics which is introduced on going from one-component to 
two-component fractional quantum Hall systems and we will have much more to say about 
them later in this Chapter. The difference in behavior between sum and difference modes is 
due to the fact that the Hamiltonian is not invariant under relative translations of the two 
layers. 

The Hamiltonian of the double-layer system in the absence of interlayer tunneling may 
be written in the following form which is convenient for calculations: 

H = nu,J2[(^liL)ai{L) + al{R)ai{R)] + 

i 

\ Y.K'PgiL)p-,{L) + V^^^p,{R)p,,{R) + 2Vi^^p,{L)p^,{R)l (37) 

where aj(L) is the Landau level lowering''^ operator for particle i in the left (L) layer, V^^ = 
yRR ig f^j^Q intra-layer Coulomb interaction, V^^ is the inter-layer Coulomb interaction, 
and Pq{X) is the density operator for layer X. [For explicit calculation we ignore the finite 
thickness of the 2D layers so that V^^ = 27reV(eg) and V^^^ = exp{-qd)Vg^^ .] In Eq. (|3) 
the Hamiltonian includes infinite constant terms corresponding to the self-interaction of each 
electron in the system. It is convenient to retain these terms so that the interaction terms 
in the Hamiltonian can be expressed in terms of density operators. Since, in all subsequent 
calculations, the Hamiltonian enters only in commutators, these non-physical constant terms 
never contribute. The operators which generate the sum and difference collective modes are 
ri'l^ = [p^°(L) ±p^°(_R)]/ \/2. With the Hamiltonian expressed in terms of density operators, 
a laborious but direct calculation makes it possible to express /^^ and s^^ in terms of s^^ 
using Eq. (|3TD . Explicit expressions for the collective mode energies and the coupling matrix 
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elements are given elsewhere. Here we only comment on some of the physically interesting 
conclusions of these calculations. 

As mentioned above the energy of the n = 1 sum mode approaches the cyclotron energy 
in the k ^ limit, in agreement with Kohn's theorem. (Since interactions are invariant 
under simultaneous translations in both layers, the proof for a single-layer system''^ trivially 
generalizes to the case of the two-layer in-phase mode.) However the energy of the n = 1 
difference mode is shifted from the cyclotron energy in the long wavelength limit where it is 
given by the following expression: 

A'-{k = 0) = hu.-l (0,gV/^/^^^(g). (38) 

Here h^^{q) is the Fourier transform of the inter-layer pair correlation: 

h'^''{q) = ^{^ML)p-,{R)\^o), (39) 

and 

/i^^(g) = ^(v&o|p.(^)P-.(^)|^o) - 1, (40) 

where N is the number of particles per layer. If the layers are uncorrelated, hLR{q) = 
hLL{(l) = and A^~{k = 0) = huc- For correlated layers hiR^q) tends to be negative, at 
least at small q since the density in the left layer will tend to be reduced when the density in 
the right layer is increased, and we can expect that A^^(A; = 0) > huOc- In fact, it is possible 
to prove that interactions (of either sign!) always increase the frequency of this mode as we 
mention below. 

By expanding the expression for /^^ at small k it can be shown that ~ k'^ whereas 
fl- = -^{k'/2)Y.qX'^h^\q) + 0{k% (41) 

The n = difference mode is dipole active. It is interesting that for the difference mode the 
interaction contributions to the dipole (oc \k\^) portions of the n = and n = 1 oscillator 
strengths are identical. By definition is positive definite so that A~(A; = 0) — huc 
must also be positive. In the single-mode approximation, the n = 1 (cyclotron) difference 
mode is always shifted to higher energy by electron-electron interactions. The situation 
is similar to that for the effect of disorder on the vibration modes of the Wigner crystal 
at strong magnetic fields where pinning of the crystal shifts both intra-Landau-level and 
inter-Landau-level modes upward by the same amount. 

To determine whether or not the n = mode is gapped, it is necessary to determine 
how behaves at small k. From general properties of wavefunctions in the lowest Landau 
level it is possible^^ to conclude that^^ ~ /c^ , whereas the behavior of the = 
difference mode structure factor depends on the difference between inter-layer and intra- 
layer correlation functions. It is possible to prove^^ that s°~ ~ fc^ at small k provided that 
intra-layer and inter-layer correlation functions separately vanish at large spatial separations. 
In general it is easy to show from the plasma arguments outlined in Section |IV| that this 
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is a property of the '^m,m,n wavefunctions. However an exception occurs for n = m. In 
this case there is no distinction between the effective plasma interactions of particles in 
the same layer and particles in different layers. The weighting of particle configurations 
depends only on the total charge densities of the two-layers and only correlations in the total 
charge densities of the two-layers go to zero at large distances. For the broken symmetry 
state '^m,m,m, it is possible^^ to show explicitly that = exp(— |Hp/2) which goes to a 
constant for A; — > 0. The SMA collective mode then goes like fc^ for this special case i/the 
ground state is approximated by '^m,m,m- (The italics above are pregnant as we will see 
later.) For the {m, m,m} states, the n = difference mode (intra-Landau level) is gapless 
at long-wavelengths. It is often the case that gapless collective modes can be identified as 
Goldstone modes associated with a broken symmetry in the ground state, and that is indeed 
the case here, as we shall see. 

We emphasize that the results discussed above follow from general sum rules and are 
independent of the approximate many-body wavefunction {^m,m',n) in terms of which we 
have framed our discussion so far. Since /^^ ~ k'^ and /°~ ~ k"^ independent of the long 
wavelength behavior of it follows quite generally that the n = sum mode has a gap and 
that the = difference mode has a gap except in the case where long-range order is present 
which results in correlation functions which do not vanish at large distances. In Fig. (|Tp we 
show results obtained for the collective mode energies of a double-layer system with a total 
Landau level filling factor = 1/2 and a layer separation d = 1.5£, close to the effective 
layer separation value for which novel double-layer fractional Hall effects have recently been 
observed. Numerical calculations^^' have established that the ground state at this value 
of d/i is accurately approximated by the {m, m, n} = {3, 3, 1} Halperin^^ wavefunction and 
we have used the correlation functions*^' ^"^ of that wavefunction to evaluate the oscillator 
strengths and structure factors. For k ^ the n = 1 sum mode (the Kohn mode) is 
unshifted by interactions while the n = 1 difference mode is shifted to higher energies as 
discussed above. The shift, which is a direct measure of inter- layer correlations, should^* be 
observable in cyclotron resonance experiments in double-layer systems. Note also that both 
sum and difference n = modes have a finite gap as expected from the above discussion. 

We are now able to compare our results for the collective modes with the Chern-Simons 
Landau-Ginzburg (CSLG) theory of the double-layer system. ^^"^^ In the CSLG random- 
phase-approximation, the sum and difference density response functions are given by: 

, Nq^/Am* , , Nq^lm^A , . 

P+{^,(1) = — — ' P-i^,(l) = — (42) 

UJ'^ — UJ% ijJ^ — ijJ_ 

where the collective mode frequencies are given by loj^ = lo^ = eB/m*c, uj^ = uJc{m — 
n)/{m + n) and m* is the effective band mass of the electrons. For the single-layer case 
the CSLG random-phase-approximation predictions are correct for the dipole active mode 
and we might have expected the same to be true in double-layer systems. The sum mode 
for double-layer systems can be clearly identified with the Kohn mode.^^ However there are 
difficulties in identifying the density difference modes in this theory. From Eq. (^) we see 
that the n = difference mode is dipole active and should have a dipole oscillator strength 
proportional to V^^h^^. One might therefore be tempted to identify uj_ with the n = 
difference mode. Then, for the case of the {m, m, m} state random-phase- approximation, 
a;_ = 0, and it would then be tempting to identify this mode with the gapless n = 
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Goldstone mode. Unfortunately, the single density- difference mode calculated within the 
double-layer CSLG theory saturates the full dipole oscillator strength Nq^/m*A. This is 
not acceptable since an excitation within the lowest Landau level can not contain explicit 
dependence on the band mass m*. The second possibility is to interpret the difference 
mode obtained in Eq. (^) as the n = 1 difference mode. In this case one is faced with the 
difficulty that the mode energy is shifted downwards from the cyclotron energy by an amount 
proportional to tUc, whereas the SMA calculations show that it should be shifted upwards 
by an amount proportional to the interlayer Coulomb energy. It has been suggested^^ that 
these difficulties can be resolved by including the mixing of the vortex excitations with the 
gaussian fluctuations in the CSLG theory. Similar difficulties arise in the fermion-Chern- 
Simons theory of the single-component z/ = 1/2 state and can be resolved in that case by 
taking sufficient care with the Landau parameters of the composite-fermion Fermi liquid. 



VI. BROKEN SYMMETRIES 

Ferromagnetic states break spin rotation symmetry since they are defined by an order 
parameter (S) which gives the magnitude and orientation of the magnetization. For the 
SU(2) symmetric case with no Zeeman term, this orientation is arbitrary. As we will see 
below, the case of a double layer system is described by a pseudospin with easy-plane 
anisotropy [U(l) or XY symmetry]. Here the magnetization vector is forced to lie in the XY 
plane in the ground state. The origin of ferromagnetism in all these systems is the Coulomb 
interaction just as it is for itinerant ferromagnets like iron. Exchange effects are particularly 
crucial in a 2DEG in a large magnetic field because the kinetic energy is quenched into 
highly degenerate Landau levels. It is advantageous to follow Hund's rule and maximize 
the spin in order to make the spatial wave function fully antisymmetric, thereby lowering 
the Coulomb energy. Since the Landau level is degenerate, this spin alignment can in some 
cases be complete since it costs no kinetic energy as it does in iron.^° 

Before considering the physical consequences of this broken symmetry, let us return to 
Table | to consider how the total spin quantum number S for a state can be determined. A 
portion of our discussion here follows that of Ref . [22] . As already mentioned, states of the 
form {m, m, m} are fully ferromagnetically aligned and have total spin S = N/2. To derive 
the spin quantum numbers for the other states in Table H we write 



S|^[Z:x]=^[$'[^]«i---a7V,Air--A 



(43) 



and use the fact that 



k,l 



(44) 



We see that 
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-T.T.<h[jmz] = sis+mz]. 
i=i j=i 



(45) 



where e{i, /c)$(- ■ ■ , Zi, - ■ ■ , Zk, ■ ■ ■) = $(■ ■ ■ , Zk, ■ ■ ■ , Zi,- ■ ■) is a label exchange operator. In 
Eq. (^^, S"^ = J^k^k ^k ^k spin raising and lowering operators, and the 

prime on the sum in Eq. (|4^) indicates that k = I is excluded. In Eq. ( ^51) S is the total-spin 
quantum number. 

To identify Halperin {m, m', n} states which are total spin eigenstates, we start with 
states where m = 1, n = and m' has any odd value. Since Y\i^j{Zi — Zj) is a Vandermonde 
determinant these orbital wave functions have the up-spin Landau level full (see table |). 
They are eigenstates of 5*^ with eigenvalue 

{N^ - Ni)/2 = Ni{m' - l)/2 = N{m' - l)/2(m' + 1) = Sm'- (46) 

Moreover S'"*'\E'i^m',o = since the up-spin Landau level is already full and there are no 
wave functions with a larger value of N-^. It follows that $i_m.o[-Z^] satisfies Eq. (^) with 
S = {N^ - Ni)/2 and hence that 

$:$:e(^, [j],)^r,m'AZ] = iVi$i,^.o[^]. (47) 
i=i j=i 

(This result can also be established by an explicit algebraic proof). Since e{i, [j])Q[Z]^[Z] = 
Q[Z]e{i, [j])^[Z] for any symmetric polynomial Q[Z], we have from Eqs. (^Tj) and ( ^5]) that 

[Z ■.X] = Sm'{S.m' + l)'^'l+2p [Z:x]. (48) 

In addition, it follows from Eq. (||) that S^^ [Z : x] = N{N/2 + 1)/2<I' [Z : x] for any 
completely antisymmetric function ^[Z], and in particular for generalized Laughlin states 
with m = m! = n. These states are merely the S"^ = members of the set of the (A^ + 1) 
fully polarized Laughlin states which are degenerate in the absence of the Zeeman term. 

We return now to the question of the physical consequences of the spontaneously broken 
symmetry of ferromagnetic states. We will focus initially on the SU(2) invariant case of 
'real' spins with zero g factor. Consider a state with v = 1 and all spins up. Because the 
up Landau level is maximally filled, the Pauli principle forces us to flip a spin if we are to 
move an electron to create a pair of charged excitations. This is illustrated in Fig. a] and 
shows that spin and charge are intimately connected in this case. 

Transport dissipation measures the thermal activation of charged excitations. In the 
absence of interactions, the energy cost of charged excitations is zero and there will be no 
u = 1 quantum Hall plateau because there is no gap. In the presence of Coulomb interactions 
a flipped spin particle-hole pair causes a loss of exchange energy of^^ 

where e is the bulk dielectric constant. This occurs because the electron spatial wave function 
is no longer perfectly antisymmetric. This energy gap is quite large (~ 150K ai B = lOT) 



17 



and is vastly larger than the bare Zeeman splitting. Hence Coulomb interactions and the 
associated ferromagnetism play a dramatic role in producing the charge gap at u = 1 (and 
because of the spontaneous magnetic ordering, will continue to do so, even if g is strictly 
zero). 

The exchange energy cost of particle-hole pair excitations is so large that it is worth 
searching for some modified form of the excitation which is less costly. A prescient analysis 
of smooth spin textures by Sondhi, et al.^^ yielded the exciting idea that 'skyrmion'^° like 
spin textures (shown in cross section in Fig. p|b]) can have relatively low energy and carry 
fermion number proportional to their topological charge (Pontryagin index)^^'^^ 

AN = -^ fdh e^, m(r) • [9^m(r) x 9,m(r)] , (50) 

where m(r) is the unit vector field representing the local spin orientation. The fermion 
number AA^ is an integer multiple of v because it is the number of times the unit sphere 
is wrapped around by the order parameter. That is, it is the winding number of the spin 
texture. For the Laughlin parent states v = 1/m, elementary spin-textures carry the same 
fractional charge as the quasiparticles discovered by Laughlin^ for spinless electrons. As we 
discuss below, the fact that the charges are the same follows from very general considerations. 
Actually, the spin texture states we have defined must contain precisely the same number of 
particles as \ipo) since the spin-rotation operator does not change the total electron number. 
However the spin-density may contain a number of well-separated textures with well-defined 
non-zero topological charge densities and hence well localized charges; only the net charge 
in the spin-texture states defined above will be zero. The system clearly has states with 
locally non-zero net charge in the spin textures. 

A simple variational wave function for a skyrmion of size A centered on the origin and 
carrying p units of topological charge is given by 

V'A = n ( ) *mmm, (51) 

where $mmm is defined in Eq.(|]), ()j refers to the spinor for the jth particle, and the 
variational parameter A is a fixed length scale. This is a skyrmion because the x — y 
component of the spin has a vortex centered on the origin and the z-component is purely 
down at the origin (where Zj = 0) and purely up at infinity (where \Zj\ ^ A) as shown 
in Fig. 0. The parameter A is simply the size scale of the skyrmion. ^^'^^ Notice that in 
the limit A — > [where the continuum effective field theory is invalid (see Section |V11| ), 
but this microscopic wave function is still sensible] we recover a fully spin polarized filled 
Landau level with p Laughlin quasiholes at the origin. Hence the number of flipped spins 
associated with the presence of the skyrmion interpolates continuously from zero to infinity 
as A increases. 

In order to analyze the skyrmion wave function in Eq. (|5T|) , we use the Laughlin plasma 
analogy.^ In this analogy the norm of tpx, Tr^a} J D[z]\'^[z]\'^ is viewed as the partition 
function of a Coulomb gas. In order to compute the density distribution we simply need to 
take a trace over the spin (we specialize here to the case p = 1 for simplicity) 
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This partition function describes the usual logarithmically interacting charge m Coulomb 
gas with uniform background charge plus a spatially varying impurity background charge 
Ap6(r), 

AP^W - -^VVW ^ (53) 



V{r) = Uog{r' + X'). (54) 

For large enough scale size X ^ i, local neutrality of the plasma^^ implies that the 
excess electron number density is precisely ^Apfe(r), so that Eq.(0) is in agreement with 
the standard result for the topological density,^^ and the skyrmion carries electron number 
1/m (for p = 1) and p/m in general. These objects are roughly analogous to the Laughlin 
quasi-hole. Explicit wave functions for the the corresponding quasi-electron objects ('anti- 
skyrmions') are more difficult to write down just as they are for the Laughlin quasi-electron 
due to the analyticity constraint.^ 

Sondhi, et al.^^ have shown that for the case of pure Coulomb interactions (i.e., with no 
finite inversion layer thickness corrections) the optimal skyrmion configuration costs precisely 
half the energy of the simple spin fiip.^^ The reason for this is simply that the skyrmion keeps 
the orientation of spins close to that of their neighbors and so loses less exchange energy. 



This will be discussed in more detail from a field theoretic point of view in Section [VII . 

Optical^^ and standard transport^^ experiments show that the charge excitation gap is 
indeed much larger than would be expected if interactions were neglected and has approx- 
imately the correct Coulomb scale, although there is not yet precise agreement between 
the observed gap and the best estimates including the finite g factor. ^^'^^ The main source 
of error is probably neglect of finite thickness effects. Calculations including finite thick- 
ness corrections do not exist at present. Quantum fluctuation effects (i.e., corrections to 
Hartree-Fock) may also be important. It should be noted that the uniform ground state 
of the V = 1 ferromagnet does not have quantum fluctuations (Hartree-Fock is exact here). 
However for flnite g factor, the length scales associated with the skyrmion are small and 
quantum corrections could well become important. 

The idea that skyrmions are the lowest energy excitations has received very strong and 
unequivocal support from numerical simulations which show that quite remarkably, adding 
a single electron to a z/ = 1 system (with g = Q) suddenly changes it from a fully aligned 
ferromagnet {S = N/2) to a spin singlet {S = 0) due to the formation of a skyrmion. ^^'^^ 
(It should be noted that this occurs in the spherical geometry. Things are slightly more 
complicated on the torus. ^^) 

The notion of charges being carried by skyrmion textures has received additional dra- 
matic experimental conflrmation in recent optically pumped NMR measurements by Barrett 
et al.^^'^'^ Their Knight shift measurements (see Fig.[^]) indicate that the electron gas spin 
polarization has a maximum at fllling factor u = 1 and falls off sharply on each side. The 
rate of fall off indicates that each charge added or removed from the Landau level turns over 
about 7 spins. This is consistent with the charge being carried by skyrmions of flnite size. 
The size of the skyrmion is determined by a competition between the Zeeman coupling which 
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wants to minimize the number of flipped spins and the Coulomb self-energy which wants to 
expand the skyrmion to spread out the excess charge density over the largest possible area. 
As mentioned above, Rezayi found numerically that a single charge entirely destroys the 
spin polarization of a system (if g = 0). Using an effective field-theory approach, Sondhi et 
al.^^ have estimated the skyrmion size in the regime of small g factor. Microscopic Hartree- 
Fock calculations expected to be more accurate in the physically accessible Zeeman energy 
regime have been performed by Fertig et al.^^ These estimates are roughly consistent with 
the experimental value of 7 spins per unit of charge. It should be noted that for finite g 
factor, the energetic advantage of the skyrmion over the simple flipped spin is considerably 
reduced. Wu and Sondhi have shown that in higher Landau levels, skyrmions have higher 
energy than other charged excitations.^^ 

For the SU(2) symmetric case discussed in this section, the existence of a vector order 
parameter (S) is in some sense trivial because the magnetization commutes with the Hamil- 
tonian. For the case of double layer systems we will see that the the pseudospin Hamiltonian 
has only U(l) symmetry and the fact that a non-zero expectation value of (S) appears is 
highly non-trivial. We will discuss this in more detail in Section |V11| . 



In this section we study the ferromagnetic broken symmetry ground state and its ex- 
citations from the point of view of (quantum) Ginzburg-Landau effective field theory. We 
will begin with the SU(2) invariant case of 'real' spin and then move on to the pseudospin 
analogy in double-layer systems. We will give here an introductory qualitative discussion of 
the physics. Part of our discussion follows fairly closely the presentation in Ref.[46]. The 
technical details of the calculations can be found there and elsewhere.^^'^^'^^ 

The standard first step in this procedure is always to identify the slowly fluctuating order 
parameter field, which in this case we have reason to believe is the local magnetization. We 
believe that on long length scales the (coarse-grained) magnetization fluctuates very slowly, 
because we know that the zero-wavevector component of the spin operator (i.e., the total 
spin) commutes with the Hamiltonian and so is a constant of the motion. We will focus on 
slow tilts of the spin orientation ignoring variations in the magnitude of the coarse-grained 
magnetization, and so define the order in terms of a local unit vector field m(r). 

General symmetry arguments can now be used to deduce the form of the Lagrangian. 
We can not have any terms that break spin rotational symmetry and thus the leading term 
which is an allowed scalar is 



where ps is a phenomenological spin stiffness coefficient and the energy is relative to the 
ground state energy. It expresses the cost due to loss of Coulomb exchange energy when the 
spin orientation varies with position. For the SU(2) invariant u = 1 case, this stiffness may 
be computed exactly. For Coulomb interactions (with no finite thickness corrections) this 
calculation yields 



VII. FIELD THEORETIC APPROACH 




(55) 




(56) 
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In two dimensions the stiffness has units of energy and is approximately 4 K at a field of 10 
T. Numerical estimates for i/ = 1/3 yield a value which is about 25 times smaller. 

As usual, there is a linear time derivative term in the Lagrangian which can be deduced 
from the fact that each spin precesses under the influence of its local exchange field. Equiv- 
alently we may note that when the orientation of a spin is moved around a closed loop, the 
quantum system picks up a Berry's phase^^ proportional to the solid angle Q enclosed by the 
path uj of the tip of the spin on the unit sphere as shown in Fig. [Q. Noting that a charged 
particle moving on the surface of a unit sphere with a magnetic monopole at the origin also 
picks up a Berry's phase proportional to the solid angle subtended by the path^*^ we may 
express the Berry's phase for a spin 5* as 

or equivalently 

gn = e^/<it%rA(m)^ (58) 

where A(m) is the vector potential of a unit monopole^'^''*^'^^ at the center of the sphere 
evaluated at the point m. That is, Vm x A = m. This phase is correctly reproduced in the 
quantum action by adding the following total derivative term to the Lagrangian for the spin 

Li = 5^-A(m). (59) 

Using the fact that the electronic density is t^/27r£^, the analog for the present problem of a 
large collection of electrons with S* = 1/2 at filling factor u may be simply written 

A = 7^/^^-^-A[m], (60) 



AttP J dt 
which yields the Lagrangian 

^ = ^/ d'r^ ■ A[m] - ip, | dV[Vm^(r)] ■ [Vm'^(r)]. (61) 

We shall see shortly that higher gradient terms can be unexpectedly significant, but this 
Lagrangian is adequate to recover the correct spin wave collective mode. Taking the spins 
to be aligned in the z direction and looking at small transverse oscillations at wave vector 
q we obtain from this Lagrangian the following equation of motion 

This yields the dispersion relation 

ncu = (63) 
which agrees with the long-wavelength limit of exact results obtained by a variety of 

99 100 

means. ' 
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At this point we have expanded the Lagrangian to lowest order in gradients and we have 
correctly found the neutral collective spin wave modes. Their dispersion is quadratic in wave 
vector just as it is for the Heisenberg ferromagnet on a lattice. However here we have an 
itinerant magnet and we have so far seen no sign of the charge degrees of freedom. It turns 
that we have to go to higher-order gradient terms in the action to see charged objects. 



We have already seen in the discussion of Fig. [g] in Section |V|, that for a filled Landau 
level, the Pauli principle forces there to be a connection between charge excitations and 
flipped spins. It turns out that the existence of a finite Hall conductivity in this itiner- 
ant magnet causes smooth spin textures to carry charge proportional to their topological 
density. One can derive this result from a Chern-Simons effective field theory, ^^'^^ or from 
microscopic considerations involving the fact that the spin density and charge density oper- 
ators do not commute when projected onto the lowest Landau level, or from macroscopic 
considerations connecting the Berry phase term to the Hall conductivity.^^' The latter is 
the least technical and the most instructive so we shall pursue it here. 

Imagine that the order parameter of the ferromagnetic system is distorted into a smooth 
texture as illustrated in Fig.^. As an electron travels around in real space along a path 
dr which is the boundary of the region F, the spin is assumed to follow the orientation of 
the local exchange field b(r) and hence traces out a path in spin space labeled u in the 
(schematic) illustration in Fig.[Q. That is, given any sufficiently smooth spin texture, we 
can write a Hartree-Fock like Hamiltonian for the electrons which will reproduce this texture 
self-consistently 

N 

H = Y.Hr).S,. (64) 

i=i 

If we drag one electron around in real space along dT and its spin follows the local b(r) 
adiabatically (which we expect since the exchange energy is so large) then the electron will 
acquire a Berry's phase Q/2 where Q is the solid angle subtended by the region u shown in 
Fig-i- 

In addition to the Berry's phase from the spin, the electron will acquire a Bohm-Aharonov 
phase from the magnetic flux enclosed in the region F. At least in the adiabatic limit, the 
electron can not distinguish the different sources of the two phases. The electron would 
acquire the same total phase in the absence of the spin texture if instead an additional 
amount of flux 

A$ = ^$0, (65) 

(where $o is the flux quantum) were added to the region F. We know however that adding 
flux to a region in a system with a finite Hall conductivity changes the total charge in that 
region. ^''^ To see this, let $(t) be the time-dependent flux inside F. Then the electric field 
along the perimeter obeys (from Faraday's Law) 

1 (i$ 

, E-rfr = --^. (66) 

dr c dt ^ ^ 

Because of the Hall conductivity (and the fact that a^x = 0), the field at the perimeter 
induces a current obeying 
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z ■ J X dr = (Txy E ■ dr. (67) 

Integrating this expression along the boundary and using the continuity equation, we have 
that the total charge inside T obeys 

-TL = — , (68) 
dt c dt' ^ ^ 

or 

Ag = -ez/— , (69) 

where we have used the fact that the Hall conductivity is quantized (and is negative for 
B = |B|z): 

o-xy = -ve^/h. (70) 

Thus V electrons flow into the region for each quantum of flux added to the region. This 
makes sense when we recall that there is one state in each Landau level per quantum of flux 
penetrating the sample. 



From Eq. (|65D we see that the spin texture thus induces an extra charge of 

Ag = -e3. (71) 

The solid angle O is of course a functional of the spin texture in the region F. For simplicity 
of analysis of this functional let us consider making up F out of a set of infinitesimal square 
loop circuits in real space of the form 

(x, y) — > (x + dx, y) — > (x + dx^ y + dy) — > (x, y + dy) — > (x, y). (72) 

The corresponding circuit in spin space illustrated in Fig.0 is, 

m(x, y) — > m(x + c/x, y) — > m(x + dx, y + dy) — > m(x, y + dy) — > m(x, y). (73) 

Approximating this circuit as a parallelogram as shown in Fig. 0, the solid angle subtended 
is (to a sufficient approximation) 

doj = [m(x + dx, y) — m(x, y)] x [m(x, y + dy) — m(x, y)] ■ m(x, y). (74) 

This may be rewritten in a suggestive form which tells us the curl of the Berry 'connection: '^^ 

du = -e^i, m ■ dj^iin x di^m dx dy. (75) 

We can now add up all the infinitesimal contributions to obtain 

Q = J^dxdy ^e^^ m ■ d^m x d^m, (76) 
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which yields a total charge of 



AQ = -— / dxdy tny m ■ d^m x di^vci, (77) 
ovr JT 

or a local charge density deviation of 

cv 

5p = -—e^^m.- d^m^i dum. (78) 
ovr 



The expression on the right hand side of Eq. ( |78D is simply the Pontryagin topological charge 
density of the spin texture. Its integral over all space is an integer and is a topologically 
invariant winding number known as the Pontryagin index. The spin textures which have non- 
zero Pontryagin index are the 'skyrmion' configurations illustrated in Fig.pb]. A microscopic 
variational wave function for these spin textures was discussed in Section |VI[ 

The charge density in Eq.(^) can be viewed as the time-like component of a conserved 
(divergenceless) topological 'three-current' which results in the following beautiful formula 

r = e"^^ eab, m''{v)dpm\v)d^m\v). (79) 

Using the fact that m is a unit vector, it is straightforward to verify that d^j^ = 0. We 
note that the fact that the expression for the topological current is not parity invariant is 
a direct reflection of the lack of parity symmetry in the presence of the external magnetic 
field. 

The mechanism we have seen here that associates charge with flux is the reason that 
quantum Hall fluids are described by Chern-Simons theories,^' '^3' ^'''^^ and is the same 
mechanism which causes Laughlin quasiparticles (which are topological vortices) to carry 
quantized fractional charge proportional to the quantized value of (Txy^ 

Having established that the electron charge density is proportional to the topological 
density of the spin order parameter field, we must now return to our Lagrangian to see what 
modifications this implies. We have already taken into account the long-range Coulomb force 
but it led only to the local spin gradient term whose coefficient is the spin wave stiffness. 
There are however additional effects of the charge (topological) density fluctuations which 
we must take into account 

/: = d'r^ ■ A[m] - ^p, | rfV[Vm^(r)] ■ [Vm'^(r)] 

q q 

Here 6pq is the Fourier transform of the charge density in Eq. ([TSD . Note that it is second 
order in spin gradients. The first of the new terms in Eq.(^) represents the coupling of 
the charge fiuctuations to the external and random disorder potentials ^(— q) and is second 
order in spin gradients. The second new term represents the mutual interaction of the charge 
fiuctuations via the Coulomb potential. Note that this is fourth-order in gradients (and so is 
not a duplicate of the ps term which also comes from the Coulomb interaction). In general 
there will be additional fourth-order terms allowed by symmetry, but we do not bother to 
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write them down since they will not have the divergent Coulomb interaction coefficient 27r/eg 
which makes the term we have kept effectively third order in q. 

We can immediately conclude several interesting things from the rather peculiar nature of 
our itinerant ferromagnet. First, unlike the case of a regular ferromagnet, a scalar potential 
can induce the formation of charged skyrmions in the ground state. Thus sufficiently strong 
disorder would have the effect of greatly reducing the net spin polarization of the ground 
state, something which should be directly observable experimentally. 

Secondly, we note that (at the classical level) the energy of a skyrmion due to the gradient 
term is scale invariant 



because we have two spatial integrations and two derivatives. [The quantity Eq is defined 
in Eq. (|49D .] Now however the Coulomb self-energy will want to expand the size of the 
skyrmion. In real spin systems this effect competes against the small but (usually) non-zero 
Zeeman coupling which wants to minimize the number of ffipped spins. This competition 
has been studied in some detail and appears to be essential to explain the experiments of 



VIII. INTERLAYER COHERENCE IN DOUBLE LAYER SYSTEMS 

The details of the double layer experiments of Murphy et al.^^ are described in the 
Chapter by J. P. Eisenstein. Here we briefly introduce the main ideas. Double layer quantum 
Hall systems (and wide single well systems^°) exhibit a variety of non-trivial collective states 
at different ffiling factors. Here we will focus on the case of total ffiling factor v = 1 (that 
is, 1/2 in each layer) which is most closely analogous to the fully ferromagnetic broken 
symmetry state for i/ = 1 with 'real' spins that we have discussed in the previous sections. 
There are many other interesting states which we do not have room to discuss here. One 
example, is the state at total filling factor z/ = 1/2 (that is, 1/4 in each layer) which is 
believed to be described by Halperin's {3,3, 1} wave function.'^' ^^'^^'^^'^^ This state is more 
nearly like a gapped spin liquid state, although, as we have already seen, it does not satisfy 
the Fock cyclic condition and so is not a true spin singlet. 

The schematic energy level diagram for the growth direction degree-of-freedom in the 
double-layer system is shown in Fig.0]. For simplicity we assume that electrons can only 
occupy the lowest electric subband in each quantum well. If the barrier between the wells is 
not too strong, tunneling from one side to the other is allowed. The lowest energy eigenstates 
split into symmetric and antisymmetric combinations separated by an energy gap Asas 
which can, depending on the sample, vary from essentially zero to hundreds of Kelvins. The 
splitting can therefore be much less than or even greater than the interlayer interaction 
energy scale, Ec = e^/{ed). 

The analogy with the spin systems studied in the previous sections is that within the 
approximations just mentioned, the electrons in a two- layer system have a double- valued 
internal quantum number — namely the layer index. The tunnel splitting in the double well 
plays the role of the Zeeman splitting for spins. 




(81) 



Barrett et al. 
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In addition to double quantum wells, there are wide single quantum wells in which the 
two lowest electric subband states are strongly mixed by Coulomb interactions. These 
systems exhibit very similar physics, and can also be approximately modeled as a double 
layer system. 

Throughout our discussion we will assume that the 'real' spins are aligned and their dy- 
namics frozen out by the small Zeeman energy. This is not necessarily a good approximation 
in experimentally relevant cases but greatly simplifies matters. Dynamics of 'real' spins in 
double layers is currently a topic of investigation. ^"^'^'^^ 



A. Experimental Indications of Interlayer Phase Coherence 

Here we will very briefiy review the main experimental indications that double well and 
wide single well systems ai u = 1 can show coherent pseudospin phase order over long length 
scales and exhibit excitations which are highly collective in nature. 

When the layers are widely separated, there will be no correlations between them and 
we expect no dissipationless quantum Hall state^''^ since each layer has z/ = 1/2. For smaller 
separations, it was predicted theoretically^^' and subsequently observed experimentally 
that there is an excitation gap and a quantized Hall plateau. The resulting phase 
diagram is shown in Fig. |p and discussed in more detail in the Chapter in this volume by J. 
P. Eisenstein. The existence of a gap has either a trivial or a highly non-trivial explanation, 
depending on the ratio Asas/-£'c- For large Asas the electrons tunnel back and forth so 
rapidly that it is as if there is only a single quantum well. The tunnel splitting Asas is then 
analogous to the electric subband splitting in a (wide) single well. All symmetric states are 
occupied and all antisymmetric states are empty and we simply have the ordinary u = 1 
integer Hall effect. Correlations are irrelevant in this limit and the excitation gap is close to 
the single-particle gap Asas (or huc, whichever is smaller). What is highly non-trivial about 
this system is the fact that the i/ = 1 quantum Hall plateau survives even when Asas -E-c 
(See Fig.0). In this limit the excitation gap has clearly changed to become highly collective 
in nature since the observed^^'^" gap can be on the scale of 10-20 K even when Asas ~ IK. 
Because of the spontaneous broken symmetry,^^"*^^' the excitation gap actually survives 
the limit Asas — ^ 0! This cross-over from single-particle to purely collective gap is quite 
analogous to the result we discussed earlier that for spin polarized single layers, the excitation 
gap survives the limit of zero Zeeman splitting. Hence, to borrow a delightful phrase from 
Sondhi et al.,^^^ 'i/ = 1 is a fraction too.' 

A second indication of the highly collective nature of the excitations can be seen in the 
Arrhenius plots showing thermally activated dissipation.^^' The low temperature activation 
energy A is, as already noted, much larger than Asas- If A were nevertheless somehow a 
single-particle gap, one would expect the Arrhenius law to be valid up to temperatures of 
order A. Instead one observes a fairly sharp leveling off in the dissipation as the temperature 
increases past values as low as ~ 0.1 A. This is consistent with the notion of a thermally 
induced collapse of the order that had been producing the collective gap. This behavior is 
very similar to that seen in 'real' spins. ^^'^^ 

The third significant feature of the experimental data pointing to a highly-ordered col- 
lective state is the strong response of the the system to relatively weak magnetic fields Bh 
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applied in the plane of the 2D electron gases. Within a model that neglects higher electric 
subbands, we can treat the electron g strictly two-dimensional. This is important 

since B\\ can then affect the system only if there are processes that carry electrons around 
closed loops containing flux. A prototypical such process is illustrated in Fig.|p. An elec- 
tron tunnels from one layer to the other at point A, and travels to point B. Then it (or 
another indistinguishable electron) tunnels back and returns to the starting point. The 
parallel field contributes to the quantum amplitude for this process (in the 2D gas limit) 
a gauge-invariant Aharonov-Bohm phase factor exp (27ri$/$o) where $ is the enclosed flux 
and $0 is the quantum of flux. 

Such loop paths evidently contribute significantly to correlations in the system since the 
activation energy gap is observed to decrease very rapidly with B\\, falling by factors of 
order two or more until a critical field, B*^ ~ 0.8T, is reached at which the gap essentially 
ceases changing. To understand how remarkably small B*^ is, consider the following. We 
can define a length L\\ from the size of the loop needed to enclose one quantum of flux: 
L||S|d = $0. (L||[A] = 4.137 X 10V(i[A]S||[T].) For B^ = 0.8T and d = 210A, = 2460A 
which is approximately six times the spacing between electrons in a given layer and more 
than twenty times larger than the quantized cyclotron orbit radius i = (hc/eB^Y^'^ within 
an individual layer. Significant drops in the excitation gap are already seen at fields of O.IT 
implying enormous phase coherent correlation lengths must exist. Again this shows the 
highly collective nature of the ordering in this system. 

The fourth indication about the nature of the coherent ordering is the fact that the 
gapped quantum Hall state at z/ = 1 can survive a finite amount of imbalance in the layer 
charge densities. Charge imbalance can be controlled by applying a gate voltage which 
increases the equilibrium charge density in one layer and decreases it in the other. In 
some cases, such as the gapped state at u = 1/2, layer imbalance immediately destroys the 
quantum Hall state. For u = 1 it does not. This is a strong hint about the different natures 
of the ordering in the two cases. We defer discussion of what this hint means to subsection 

iviTra. 



B. Effective action for double layer systems 

Having established the correct form of the effective low-energy Lagrangian for the SU(2) 
invariant case, let us now turn to the U(l) symmetric case in double layer systems. The spin 
analogy is relatively straightforward, but can be confusing until one gets used to it. The 
idea is simply that each electron can be in either the upper layer or the lower layer and we 
refer to these two states as pseudospin up and down respectively. For the moment we shall 
assume the layers are identical and also neglect the possibility of tunneling between the two 
layers. The confusing point is that quantum mechanics nevertheless forces us to consider 
the possibility that an electron can be in a coherent superposition of the two pseudospin 
states so that its layer index is uncertain. 

Portions of our discussion here and in the following subsections follows that of Ref.[46]. 
Further technical details can be found therein. The formal mapping that we use to define the 
pseudospin density operators is the following. The z component of the pseudospin density 
represents the local charge density difference between the layers 
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(82) 



The X and y components of the pseudospin density are off-diagonal and can be combined to 
form tunnehng operators 



S+{r) = S-{r) 



1 



^}(r)^^(r). 



If for instance, tunnehng is present then the Hamiltonian contains a term 



T^-tJ d^r [S+{v) + S-{r)\ = -2t J dV5^(r). 



(83) 



(84) 



We know that ferromagnetism in real spin systems is not the result of truly spin- 
dependent forces but rather a byproduct of Coulomb exchange forces. Hence the effective 
spin Hamiltonian must contain only spin scalars such as Sj ■ Sj. Here however the Coulomb 
forces are explicitly pseudospin dependent since the intra- and inter-layer Coulomb interac- 
tions are not identical. Following our previous discussion^^ we define 
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V.'^liV.'^-V,^), (86) 

where V,f is the Fourier transform with respect to the planar coordinate of the interaction 
potential between a pair of electrons in the same layer and Vjf is the Fourier transform of 
the interaction potential between a pair of electrons in opposite layers. If we neglect the 
finite thickness^^ of the layers, = 27re^/(eg) and = cxp{—qd)V^, where d is the 
layer separation. The interaction Hamiltonian can then be separated into a pseudospin- 
independent part with interaction and a pseudospin- dependent part. The pseudospin 
dependent term in the Hamiltonian is 

k 

Here Sk is the Fourier transform of the spin density at wave vector k and the overbar 
indicates projection onto the lowest Landau level. Since V/^ > Vjf, is positive and 
this term produces an easy-plane, as opposed to Ising, pseudospin anisotropy. That is, this 
term prefers for the spin to lie in the XY plane. If the spin orientation moves up out of 
the XY plane so that (S*^) 7^ 0, then the energy increases. We can view this energy cost as 
simply the charging energy of the capacitor formed by the two layers since the pseudospin 
component measures the charge difference between the two layers. 

The pseudospin symmetry of the Hamiltonian is reduced from SU{2) to U{1) by this 
term. In addition, this term increases the quantum fiuctuations in the system since it does 
not commute with the order parameter 

[Fsb,^l7^0, (88) 
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where = x,y. Thus total spin is no longer a sharp qTiantum number. However for 
small layer separations, we expect the quantum fluctuations to be small. At very large 
layer separations we expect the quantum fluctuations to become dominant and produce the 
disordering phase transition which uncouples the two layers. 

In the absence of the symmetry breaking term {d — 0), we know the exact quantum 
ground state for z/ = 1 since it is simply a fully occupied pseudospin-polarized Landau level. 
This state is 25* + 1 = + 1 fold degenerate (since there is no charging energy in this 
unphysical limit). As a first approximation for finite d, we assume that the eigenstates 
remain exactly the same, but the charging energy lifts the degeneracy among these states by 
favoring the ones with 5"^ ~ 0. That is, we assume that the spin vector is still fully polarized 
and non-fluctuating, but it now lies in the XY plane. 

We argue that the form of the energy-functional we derive must remain valid even when 
quantum fluctuations due to the pscudospin-dependent terms in the Hamiltonian are present. 
However, the coefficients which appear in the energy-functional will be altered by quan- 
tum fluctuations and the explicit expressions we derive below are accurate only when the 
pseudospin-dependcnt interactions are weak, i.e., only when the layers are close together. 
Estimates have been obtained for the quantum fluctuation corrections to these coefficients 
from finite-size exact diagonalization and many-body perturbation theory calculations.^^' 

To better understand what we are assuming, consider the specific example of a state 
which has S — N/2 and has its spin oriented in the x direction so that it is an eigenstate of 



Here X is a LLL orbital label and |0) is the electron vacuum. It is clear from the construction 
that this state has every orbital in the LLL filled so it has v = 1. However each electron 
is in a linear combination of up and down pseudospin which puts it in an eigenstate of 
S^. This particular combination is familiar because it is the symmetric state which is the 
exact (non-interacting) ground state in the presence of tunneling between the layers. In the 
absence of V^b this would be an exact eigenstate of the system (with or without tunneling) 
even in the presence of the symmetric interaction term Vq because it is a filled Landau level. 
This state has zero net charge on the capacitor on the average {{S^) = 0) and so first-order 
perturbation theory in Vsb favors it for the ground state. 

The confusion that many people have at this point is over the fact that this is a coherent 
state with uncertain S^. It makes perfect sense if the Hamiltonian includes tunneling that 
the charge in each layer will be uncertain. However we are arguing that this state is a good 
approximation to the ground state even in the absence of tunneling because it has favorable 
Coulomb energy. That is, we are assuming that the system spontaneously breaks the U{1) 
symmetry associated with conservation of layer charge difference^^"^^' and is acting much 
like a superfiuid^^ which breaks U{1) symmetry. We refer to this state as having spontaneous 
interlayer phase coherence. It is quite analogous to a BCS state which has an uncertain total 
number of Cooper pairs. In the absence of tunneling, the electron energy cannot be sensitive 
to the actual relative phase between the two layers and so any coherent state of the form 



total S"" 
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will have equally good energy. The global U(l) phase simply determines the orientation 
of the total spin in the XY plane 

(^I^+I^) = (^1^^ + iSy]"^) = —e'^. (91) 

This is exactly analogous to the situation in a BCS superconductor where the energy is 
independent of the phase which determines the coherence between states of different 
numbers of Cooper pairs 

(^{^{) ~ e^. (92) 

The complex order parameter in a double layer system given by Eq.(0) has an amplitude 
and a phase like that of a superconductor, but in many ways is more reminiscent of an 
excitonic insulator since it is charge neutral (i.e., it contains not '?/'j'?/'|) and physically 
represents a particle in one layer bound to a hole in the other layer (but we do not know 
which layer contains the particle and which contains the hole). Some years ago, Datta^^° 
considered similar states in double layer systems in zero B field. 

At first, it may seem counter-intuitive that, in the absence of tunneling, the system 
energy could depend on the relative phase for finding the particle in the upper or lower 
layer. Indeed the energy is unchanged when a constant is added to However the energy 
does depend on gradients of The 'spin stiffness' ps is non-zero because of the loss of 
exchange energy which occurs when varies with position. Imagine that two particles 
approach each other. They are in a linear superposition of states in each of the layers (even 
though there is no tunneling!). If these superpositions are characterized by the same phase, 
then the wave function is symmetric under pseudospin exchange and so the spatial wave 
function is antisymmetric and must vanish as the particles approach each other. This lowers 
the Coulomb energy. If a phase gradient exists, then there is a greater amplitude for the 
particles to be near each other and hence the energy is higher. This loss of exchange energy 
is the source of the finite spin stiffness and is what causes the system to spontaneously 
'magnetize'. We may describe this by a term in the effective action of the form 

E = \p,\d\\^^\\ (93) 

We see immediately that at finite temperatures this system will be described by a classical 
XY model and will undergo a Kosterlitz-Thouless phase transition. We will cover this in 
more detail further below. 

C. Superfluid Dynamics 

We consider now the dynamics of a two layer system with easy-plane pseudospin 
anisotropy but no tunneling. Let us back up and note that implicit in Eq.(P3D is the 
assumption that the spin lies exactly in the XY plane. To be more complete we should 
derive the effective Lagrangian in the presence of the easy plane anisotropy. On general 
symmetry grounds this must be of the form (neglecting higher order gradients needed to 
describe charge fluctuations) 
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£ = I dPr^ . A[m] - j Sr + ^|Vm^|^ + ^[|Vm^|2 + |Vm^p]} . (94) 

Within the Hartree-Fock approximation which is presumably vahd only for small layer sep- 
arations, the coefficients P,Pa, and pe can be evaluated in terms of the density- density 
correlation function of the ground state. (The gradient expansion in is not strictly 
valid since it turns out^^ that the long range of the Coulomb interaction leads to a non-local 
m^m^ interaction term not included above. This term vanishes in the SU(2) invariant case 
but here is more important than the |Vm^p term and less important than the /5|m^p term 
at long wavelengths. We retain the |Vm^p term in the following expressions only to remind 
us of its importance in the SU(2) invariant limit.) 

Since is massive, the equations of motion derived from this Lagrangian lead to a 
linear rather than quadratic collective mode dispersion^^"^^'^^'^^''^"^ like that of the Goldstone 
mode in a superfluid or an antiferromagnet. We take the pseudospin of the system to be 
polarized in the x direction and consider the linear response to a time- and space-dependent 
pseudospin 'Zeeman' field in the y — z plane. Using the equations of motion determined 
from the Lagrangian in Eq. (|94D and Fourier transforming with respect to both time and 
space we find that 



-too -^{2(3 + q^pA) 
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where hy and hz are the Fourier coefficients of the pseudospin magnetic field at frequency 00 
and wavevector q. Physically hz corresponds to a time- and space-dependent bias potential 
between the two wells, while hy could arise from a space- and time-dependent interlayer 
tunneling amplitude. We see immediately that the response has a singularity at the collective 
mode frequency 
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+ q'ps. (96) 



For the djt = case, j3 = 0, pA = Ps = Ps, and the collective mode frequency reduces to 
the result obtained previously for the spin-wave collective mode of isotropic ferromagnets 
47rg^p°/i/]. The collective mode corresponds to a spin-precession whose ellipticity 



UJ. 
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increases as the long-wavelength limit is approached. The presence of the mass term {f3 7^ 0) 
changes the collective mode dispersion at long wavelengths from quadratic to linear. In 
the limit of small q 

4:71 



uj, = —p(3psq. (97) 

Thus we see that the system is acting very much like a bosonic superfluid of weakly 
repulsively interacting particles which has a Goldstone mode (if and only if there is non-zero 
repulsion). ^^"^^'^^'^^ However we must again emphasize that we have a charge neutral order 
parameter. The singular response occurs for electric fields of opposite signs in the two layers. 
Another way to say this is that the 'charge' conjugate to the U(l) phase field is . The 
current associated with phase gradients. 
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Jzz = ^V^, (98) 

is the flow of pseudospin density wliicli is tlie difference between the electrical currents in 
the two layers.^^'"^^'^^'^^ This simply reflects the fact that the object which condenses is a 
neutral excitonic composite of a particle and hole. 

These analogies with superfluidity naturally raise the question of analogs of other phe- 
nomenon known to occur in superfluids and superconductors. It is believed that there is 
no analog of the Meissner effect in these systems. ^'^"^^ However, it has been suggested that 
Josephson or 'Josephson like' effects should occur in the tunneling transport between the 
two layers. ■^^"'^^ In this picture one views each layer as analogous to one side of a Josephson 
junction. This appealing analogy appears to imply measureable consequences only in the 
case where superconducting contacts are made to each two-dimensional electron layer; the 
distinction between the system and an ordinary Josephson junction system then seems to 
be artificial. More interesting, in our view, is the recent 'sideways tunneling' proposal of 
Wen and Zee^^^ which appears to provide a quite precise analog of the Josephson effect, but 
in this case for pseudospin superfluidity. Here one views the pair of layers as constituting a 
single superfluid on one side of a Josephson junction. A separate pair of layers is imagined to 
be coupled to the first by weak 'sideways' tunneling. In each layer, there is a phase variable 



ip and a pseduospin supercurrent given by Eq.(98). When this current reaches the junction. 



it can coherently jump across the weak link and continue onwards. Physically this coherent 
tunneling of pseudospin involves an electron tunneling one way between the upper layers 
and another electron tunneling in the opposite direction between the lower layers. It should 
be possible in principle to see AC Josephson oscillations of the current at a frequency of 
2eV/h where V an appropriately applied bias voltage (of opposite sign in the two layers). 
It is unlikely however that a sample with no tunneling within each pair of layers, but finite 
tunneling across the junction, can ever be produced. One could imagine using a parallel B 
field to effectively shut off the tunneling within each pair, but even weak disorder will pin 
the discommensurations and probably ruin the effect. 

As an aside to the question of spin channel superfluidity, we point out that double layer 
systems are ideal for probing electron-electron interactions via mutual drag effects in trans- 
port at zero 

field. More details on some of the mutual transport effects associated with 
spontaneous interlayer phase coherence are presented in Ref.[46]. In addition, it has been 
predicted that for double-layer fractional quantum Hall states without broken sjTiimetries, 
but with inter-layer correlations of the type described by the Halperin {m, m', n} wavefunc- 
tions, anomalously large (and quantized!) mutual drag effects will occur. ^^'^^ 

We are now in a position to return to a discussion of the fourth experimental indication 



of interlayer phase coherence that was briefly alluded to at the end of subsection |VIII A| , 
namely the fact that unlike the z/ = 1/2 state, the v = 1 state survives a finite amount of 
charge imbalance in the two layers. The rigidity of the z/ = 1/2 state can be understood 
macroscopically from the Chern Simons effective field theory point of view,^^'^^^ or micro- 
scopically from the fact that the Halperin {3, 3, 1} state is believed to give a good description 
of the state. The plasma analogy shows that this state has a sharply defined total filling 



z/ = 1/2 and that the density in each layer must be exactly equal (see Section [TVl). Another 



way to say this is that the {3, 3, 1} state is almost a singlet spin liquid which must have 
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= (although strictly speaking, it is not an eigenstate of total spin). That is to say, it 
is pseudospin incompressible. 

In contrast to this, the u = 1 state is pseudospin compressible. This is clear from the 
existence of pseudospin superfluidity with its gapless collective Goldstone mode. The state 
is ferromagnetically ordered and there is a whole family of states differing only in 5*^ which 
would be degenerate were it not for the charging energy. The charging energy lifts this 
degeneracy and picks out the 5*^ = state. Application of a bias voltage to unbalance the 
layers, simply picks a different state having non-zero out of the manifold. 

We can supplement this picture by thinking macroscopically of the XY model describing 
the ordering. Layer imbalance corresponds to tilting the order parameter slightly up out 
of the XY plane. This changes no symmetries in the problem, and simply renormalizes the 
spin stiffness slightly because the projection of the spin onto the XY plane is reduced. Thus 
the ordering which produces the charge gap is weakened but not (immediately) destroyed 
by layer imbalance. One can also describe this microscopically by modifying the variational 
wave function in Eq. (|90D to the form 

l^) = n;7^(4T+7e^^cy |0), (99) 

which obeys 

(vI/|S1vI;) = ^(l-f). (100) 

We close by noting that a class of experiments that would be very useful in probing the 
pseudospin superfluidity would be inelastic (Raman) light scattering to try to detect the 
gapless collective excitation mode. To see this at finite frequency requires finite wavevectors 
which in turn would require some sort of grating coupler. These are difficult to fabricate 
on sufficiently short length scales. Here however, it turns out that it is possible to have 
the system act as its own grating by tilting the applied magnetic field which will induce a 
twisting of the order parameter^^' as we will discuss in Section 0. 



D. Merons: Charged Vortex Excitations 

Continuing the superfiuid analogy, we study vortex excitations in this section and dis- 
cuss the Kosterhtz-Thouless phase transition induced by the unbinding of these topological 
defects.^^'^^'^^'"'^ The order parameter of the system in the presence of a vortex at the origin 
has the following approximate form 

m = I ±^Jl - [m^(r)]2cos6', - [m^(r)]2 sin 6* , m^(r)| . (101) 

Here the ± refers to right and left handed vortices respectively, and 6 is the azimuthal angle 
made by the position vector r. At asymptotically large radii, vanishes to minimize the 
charging energy. However in the vortex core we must have — > ±1 (and m^,my — > 0) 
to prevent a singularity in the gradient energy. Thus there are four fiavors of topologically 
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stable objects which we refer to as 'merons' since it turns out that they are essentially half 



skyrmions. These are illustrated in Fig. |T0 



The local topological charge density calculated from 5p = — g^ejj(9jm x djin) ■ m can be 
expressed in the form 

and the total charge is 

Q = J d\6p{r) = ^[m'{oo)-m'{0)]. (103) 

For a meron, the spin points up or down at the core center and tilts away from the z direction 
as the distance from the core center increases. At asymptotically large distances from the 
origin, the spins point purely radially in the x — y plane. Thus the topological charge is ±^ 
depending on the polarity of core spin. The general result for the topological charge of the 
four meron flavors may be summarized by the following formula: 

g = ^K(oo) -m^(0)]nv, (104) 

where is the vortex winding number. The formulae derived above for the meron charge 
do not rely for their validity on the variational ansatz assumed in Eq.( |101| ). They are quite 
general and follow from the fact that a meron topologically has half the spin winding of a 
skyrmion. The meron charge of ±1/2 is a topological invariant and implies that the electrical 
charge is ±z/e/2. 

The fact that merons carry fractional charge ±z/e/2 can be deduced from a Berry phase 
argument similar to the one used to find the skyrmion charge. We simply note that an 
electron moving at a large distance around a meron will have its spin rotated through 2tt 
in the x — y plane due to the vorticity. We know that the Berry's phase for rotating a 
spin one-half object in such a way is exp {i27cS) = —1. Thus the meron produces the same 
Berry's phase as half a flux quantum. From Eq.(|5^ we then obtain AQ = ±i>e/2. [The 
ambiguity of the sign of the charge associated with half a flux quantum can only be resolved 
by examining the behavior of in the meron core. It depends on whether the mid-gap 
state induced by the topological defect (discussed below) is empty or occupied.] 

It is instructive to write down explicit microscopic variational wave functions for vortices 
(merons). We start with the simplest example: a meron with vorticity +1 and charge — | 
that has the smallest possible core size: 

= (^y^^T + |0). (105) 

Here |0) is the fermion vacuum, c|„|(|) creates an electron in the upper (lower) layer in the 
angular momentum m state in the LLL, and M is the angular momentum quantum number 
corresponding to the edge. The vorticity is +1 because far away the spin wave function is 
essentially 
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X{0) = ^ \ , (106) 



where 9 is the polar angle. The charge is — | because we have created a hole in the center 
of the lower layer (every state has occupancy 1/2 except m = | which is unoccupied). 
Since the spin is pointing up at the center, this agrees with the spin-charge relation derived 
earlier. From the spin-charge relation we know we can change the sign of the charge of a 
meron by changing the direction of spins in the core region without changing the vorticity. 
This can be seen explicitly from the wave function: 

M 1 

= cl n (^4t + y|4+u)|0). (107) 

This state has charge +| because we have put in an electron in the m = state in the lower 
layer. Obviously what we did (in terms of the spin texture language) is to flip the spins in 
the core region to the down direction without changing the vorticity of the meron at long 
distances. In this construction, one sees that in a sense, the merons are like fractionally 
charged 'mid-gap states' in polyacetylene which can be empty or occupied. 

A meron with vorticity -1 is readily obtained by simply interchanging the labels m 
and m + 1 in the subscripts in Eq. (|105|) and Eg. (|107|) . Invariance under pseudospin reversal 
guarantees the equality of the energies of corresponding ± vorticity merons. However the two 
charge states for a given vorticity are not necessarily degenerate, just as Laughlin quasiholes 
and quasiparticles are in general non-degenerate. [For the special case u = 1 particle-hole 
symmetry guarantees degeneracy.] 

It is instructive to attempt to find a variational wave function for a pair of merons to 
show that a meron is half a skyrmion. Consider the situation where we have a pair of merons 
of opposite vorticity but the same charge, located at points Zi and Z2. To achieve this one 
meron must have = +1 in its core and the other must have = — 1 in its core. The 
following wave function seems to do the job: 





where ^mmm is defined in Eq.(^, ip is an arbitrary constant phase angle, and ()j refers to 
the spinor for the jth particle. Note that at very large distances away from Zi and Z2 the 
spinor for each particle asymptotically approaches 

(109) 

and so corresponds to a fixed spin orientation in the xy plane at an angle from the x 
axis. Hence the net vortex content of the pair of excitations is zero. It is clear from the 
construction that the spin orientation is purely up for an electron located at Z2 and purely 
down at Zi. Furthermore the net charge must be +z/e since asymptotically the factor of Zj 
is the same one as for the Laughlin quasihole in a spin polarized state. By symmetry, it 
seems that there must be half that amount, +1^6/2 associated with each of the two objects 
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located at Z\ and Zi- (It turns out however, that the charge is not spht in half and localized 
near and Z2 as we would really want for a meron pair.) 

We can easily show that this particular wave function actually is just another way to 
represent a skyrmion. Let Zx = \ and Z2 = —A and, for simplicity, assume that the spins 
are asymptotically oriented in the x direction so that = 0. Now perform a global rotation 
of all the spins about the y axis by an angle of —it/ 2. Using 

we see that we have recovered the variational skyrmion wave function studied previously in 
Eg. (pTD ■ The form in Eq.( |108| ) is the appropriate one in the U(l) symmetric case since it 



keeps the spins primarily in the XY plane. We expect that a pair of merons can be deformed 
into a skyrmion and each meron can be properly viewed as carrying half the topological and 
electrical charge of a skyrmion. However there must exist variational wave functions which 
are better for the U(l) case than the one discussed here in the sense that the present one 
has a smooth charge distribution centered in between Zi and Z2 rather than being more 
closely associated with the two points defining the meron positions. 

Finally we note that merons carry fractional statistics 1/4; they are 'quarterons'.^^ This 
can be seen from the Chern Simons theory^^ or from the fact that two of them together 
make a skyrmion which is a fermion. 



E. Kosterlitz-Thouless Phase Transition 

It is well known that the presence of vortex topological defects can disorder the XY phase 
of the ground state. This will occur even at zero temperature due to quantum fluctuations 
if the layer separation exceeds a critical value d*. One could also conveniently tune through 
this transition in a single sample with fixed d by varying a gate voltage to induce layer 
imbalance which will renormalize the spin stiffness. We focus here however on the effect of 
finite temperature and thermally induced vortices in samples which are otherwise ordered 
at zero temperature (because d < d*). 

Integrating out the massive fluctuations and going to finite temperature, we are led 
to a classical XY model of the form shown in Eq. (p3D . Hartree-Fock estimates of the spin 
stiffness at finite layer separation give values in the range 0.1-0.5 K for typical experimental 
sample parameters. The Kosterlitz-Thouless phase transition is an unbinding of vortex 
defects in the XY model and occurs at a temperature given approximately by the value of 
the bare spin stiffness. These unbound vortices will cause the long wavelength effective spin 
stiffness in Eq.(^) to renormalize discontinuously to zero at the Kosterlitz-Thouless critical 
point. 

In the usual way, the scale-invariant classical action of Eq.(^) leads to a logarithmic 
interaction among vortices. A gas consisting of M merons will have an energy of the form 

E = ME^ore - ninflirps In ( — j + ^ 1i^iY~R~' ^^^^^ 
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Here i^core is the meron core energy/^^ -Rcore is the meron core size, Rij is the separation 
between the ith and jth merons, is the vortex charge (winding number, ±1) of the ith 
meron. The last term is new and is unique to the present problem. It represents the Coulomb 
interaction among the fractional charges bound to the merons. qi = ±1 is the sign of the 
electrical charge (±e/2) of the zth meron. 

The origin of the logarithmic interaction in a superconducting or superfluid film is the 
kinetic energy stored in the supercurrents circulating around the vortices due to order pa- 
rameter phase gradients. Here the logarithmic interaction arises from the loss of Coulomb 
exchange potential energy in the presence of phase gradients produced by the vortices. 

The Coulomb interaction in Eg. ( |1 1 1| ) falls off more rapidly than the log interaction (it is 



effectively one order higher in derivatives than the log interaction) and so is perturbatively 
irrelevant at the Kosterlitz-Thouless critical point. That is, the Kosterlitz-Thouless temper- 
ature may be shifted somewhat but the phase transition itself is unaffected. However in the 
limit of strong Coulomb effects (due to small e or large vortex fugacity) the global phase 
diagram becomes extremely rich. Among other things there is a phase transition to a chiral 
state with non-zero order parameter (uiqi) in which vortex charge and electrical charge are 
no longer independent. The rich physics and the novel phase diagram of this model has been 
elucidated in an interesting paper by Tupitsyn, Wallin, and Rosengren.^^^ 

If found experimentally, this KT transition would be the first finite temperature phase 
transition in a quantum Hall system. All other transitions between plateaus, etc. are zero 
temperature transitions because the vortices there (Laughlin quasiparticles) do not interact 
logarithmically and are unconfined by an analog of the Anderson-Higgs mechanism.^' ^'^^'^^ 
As noted above the characteristic energy scale for ps is 0.1-0.5 K. This is not a problem from 
the point of view of experimentally achievable temperatures. However, in order to have the 
needed U(l) symmetry it is essential (as we discuss in Section pX[ ) that the tunneling am- 
plitude between the layers be much smaller than this scale. Such samples have not yet been 
constructed, but could be in principle because the tunneling amplitude falls off exponentially 
with layer separation while the Coulomb interactions which control the stiffness ps fall off 
only as a power law. Nevertheless it will not be easy since present samples are already close 
to the critical separation d*/i. The latter could perhaps be circumvented by going to lower 
density samples and correspondingly lower magnetic fields. 

The standard experimental signatures of the KT transition are zero linear response dis- 
sipation below Tkt, but zero critical current, and a characteristic jump in the exponent 
associated with the non-linear response 

V ~ P, (112) 

from p = 1 above the transition to p = 3 just below the transition. Dissipation is caused 
by phase slips associated with the motion of unconfined vortices. A transport current ex- 
erts a magnus force whose direction depends on the sign of the vorticity and which moves 
the vortices at right angles to the current. Below the transition the vortices are confined 
into (vorticity) neutral pairs by their logarithmic attraction and no longer couple to the 
supercurrent. 

In a superconductor or superfluid at finite temperature there is always some sort of 
'normal' fluid present (thermally excited Bogoljubov quasiparticles, for example) which can 
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produce dissipation. However there are no electric fields to couple to these particles since 
the supercurrent shorts them out. Thus there is no dissipation. Analogous normal fiuid 
excitations exist in the quantum Hall systems we are studying here. If we define normal 
fiuid as anything which is charged but has no net vorticity then we see that a pair of opposite 
vorticity, but like-charged merons constitutes normal fiuid that can be thermally activated. 
The validity of the meron pair picture has been confirmed by numerical calculations. In the 
pseudospin superfiuid channel (i.e., opposite electrical currents in each layer) this normal 
fiuid will not couple to the supercurrent since it is vortex neutral. The linear response 
dissipation will therefore drop to zero below Tkt and all the other usual signatures of the 
Kosterlitz-Thouless transition will be present. 

Doing transport experiments in the pseudospin channel requires separately contacting 
the two layers and having extremely high tunneling resistance between the layers without 
moving them so far apart that they decouple. It would be easier to do the experiment in the 
ordinary charge channel in which both layers are contacted simultaneously and current fiows 
in the same direction in both layers. Unfortunately in this channel there is a net electric 
field induced by the u = 1 quantized Hall resistivity and the normal fiuid will couple to this 
producing thermally activated dissipation.^^ While we would expect to see the dissipation 
decrease rapidly towards zero below Tkt, it would not show any discontinuity as it would 
in the pseudospin channel. This is the same mechanism that causes the small but non-zero 
thermally activated dissipation in ordinary quantum Hall plateaus. 

The energy for an excitation of 'normal fiuid' consisting of a pair of like-charged but 
opposite vorticity merons given by Eg . (|1 1 1|) is minimized at a separation of 



e2 



Ro = - . (113) 

For typical values of ps, Rq > lOi, and so is large enough to justify the field-theoretic 
continuum approximations used in deriving Eq. (|lllD . 



IX. TUNNELING BETWEEN THE LAYERS 

A finite tunneling amplitude t between the layers breaks the U{1) symmetry 



(114) 



by giving a preference to symmetric tunneling states. This can be seen from the tunneling 
Hamiltonian 

H^ = -tj d\ {V^}(r)z^i(r) + ^I(r)^t(r)} , (115) 
which can be written in the pseudospin representation as 

= -2t f SrS'^ir). (116) 



(Recall that the eigenstates of S^' are symmetric and antisymmetric combinations of up and 
down.) 
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We can shed some more light on the spontaneous symmetry breaking by considering 
the tunnehng Hamiltonian i^x in Eg .( 11151 ) as a weak perturbation. Naively, since particle 
number is separately conserved in each layer for t = 0, one might expect 

\nn^^{^\H^\ij) = 0. (117) 

That is, one might expect that the first-order term in the perturbation series for the energy 
shift due to t to vanish. Instead however, we find that the energy shifts linearly in t 



lim lim — - (-ih \Hnr\ -ijj) = lim lim 7 (i' 



= -m\ (118) 

where A is the system area, and rri^ is, by definition, the magnetization which is the system's 
order parameter. If the interlayer spacing d is taken to be zero, one can readily show^^ that 

the variational wave function in Eq.(|5UD is exact, hence lim^ ,Qm^ = 1, and t = Asas/"^- 

For finite d, Eq. (^) is no longer exact and quantum fluctuations will^^ reduce the magnitude 
of and we must renormalize the hopping parameter t appropriately. 

As the layer separation d increases, a critical point d* will be reached at which the 
magnetization vanishes and the ordered phase is destroyed by quantum fluctuations.*^^' ^^'^^ 
This is illustrated in Fig. [|| and discussed in more detail from the experimental point of 
view in the Chapter by J. P. Eisenstein.^^ For finite tunneling t, the collective mode becomes 
massive and quantum fluctuations will be less severe.^'' Hence the phase boundary in Fig.]^ 
curves upward with increasing Asas- For As as = the destruction of long-range order and 
the charge excitation gap are intimately related and occur simultaneously at d* and zero 
temperature. For finite A^^g the system always has non-zero even in the phase with 
zero charge gap. 

We have already seen that finite layer separation reduces the pseudospin symmetry from 
SU(2) to U(l). The introduction of finite tunneling amplitude destroys the U(l) symmetry 
and makes the simple vortex-pair configuration extremely expensive, thereby destroying the 
KT transition. To lower the energy, the system distorts the spin deviations into a domain 
wall or 'string' connecting the vortex cores as shown in Fig.|ni]. The spins are oriented 



in the x direction everywhere except in the shaded domain line region where they tumble 
rapidly through 2it. The domain line has a fixed energy per unit length and so the vortices 
are now confined by a linear 'string tension' rather than logarithmically. We can estimate 
the string tension by examining the energy of a domain line of infinite length. The optimal 
form for a domain line lying along the y axis is given by 

V?(r) = 2 arcsin [tanh (x/^)], (119) 

where the characteristic width of the string is 



The resulting string tension is^^^ 



(120) 
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Provided the string is long enough (i? ^ the total energy of a segment of length R will 
be well-approximated by the expression 

E'^^^ = 2E;, + ^ + ToR. (122) 

The prime on Emc in Eq.( p.22| ) indicates that the meron core energy can depend /S.sas- -^pair 

is minimized aX R = R'q = \J / (4eTo). Note that apart from the core energies, the charge 

gap at fixed layer separation (and hence fixed ps) is proportional to Tq^^ oc t^/'^ ~ ^^sas 
which contrasts with the case of free electrons, for which the charge gap is proportional to 
^SAs- Note that because the exponent 1/4 is so small, there is an extremely rapid initial 
increase in the charge gap when tunneling first becomes important. (See below.) Note that 
the present classical analysis ignores the stabilizing effect of the tunneling- induced gap on 
the quantum fluctuations of the pseudospin magnetization. 

The crossover between the meron-pair pseudospin texture which holds for t = and the 
domain line string pseudospin texture described above occurs at a finite value of t which 
we can estimate by the following argument. Recall from Eq. ( P-13| ), that the equilibrium 
separation of a meron pair in a system with no tunneling is defined to be Rq. For R'q > Rq 
the vortices are already bound by the logarithmic attraction due to the gradient energy 
before the linear attraction due to the hopping becomes important at larger separations. In 
this regime tunneling does not play an important role in determining the nature of the lowest 
energy charged pseudospin texture. As t increases Rq oc t~^^^ decreases and will eventually 
reach Rq which is, of course, independent of t. Since 



R'q _( 27rV. V^'_ 



Ro W/ieOj 4i?; 



(123) 



the characteristic width of the domain line becomes comparable to Rq in the same range 
of t values where Rq and Rq become comparable. We may conclude that the nature of 
the charged pseudospin texture crosses over directly from the meron pair form to the finite 
length domain line string form for ps/[e^/(e^)] ~ 1/25, or equivalently for t ~ tcr where 



3.9 X 10^ 



(eVe£) 



1 3 



eV(e£) 



(124) 



The crossover tunneling amplitude is thus typically smaller than 5 x 10^^[e^/(e£)]. Typical 
tunneling amplitudes in double-layer systems are smaller than ~ 10~^[e^/(e£)] and can 
be made quite small by adjusting the barrier material or or making the barrier wider. 
Nevertheless, it seems likely that t will be larger than tcr except for samples which are 
carefully prepared to make t as small as possible. As t increases beyond tcr, R'q will continue 
to decrease. When R'q becomes comparable to the microscopic length, the description given 
here will become invalid and the lowest energy charged excitations will have single-particle 
character. However, the domain-wall string picture of the charged pseudospin texture has a 
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very large range of validity since R'q oc t decreases very slowly with increasing t at large 
t. Writing R'q ~ [eV(e87rp,)](tcr/t)^/^ we find that ~ ^ only for t ~ 10-2[(eVe£)Vp,]. 
Using typical values of ps we see that the charged excitation crosses over to single particle 
character only when the hopping energy t becomes comparable to the microscopic interaction 
energy scale. The various regimes for the charge excitations of double-layer systems are 
summarized in Table fT^. Almost all typical double-layer systems lie within the regime of 
the domain-wall-string pseudospin texture charge excitation, and hence will not exhibit a 
true Kosterlitz-Thouless phase transition. 

X. PARALLEL MAGNETIC FIELD IN DOUBLE LAYER SYSTEMS 

Murphy et al.^^ have shown that the charge gap in double layer systems is remarkably 
sensitive to the application of relatively weak magnetic fields B^\, oriented in the plane of 
the 2D electron gas. Experimentally this field component is generated by slightly tilting 
the sample relative to the magnetic field orientation. Tilting the field (or sample) has 
traditionally been an effective method for identifying effects due to (real) spins because 
orbital motion in a single-layer 2DEG system is primarily^^ sensitive to B±, while the (real) 
spin Zeeman splitting is proportional to the full magnitude of B. Adding a parallel field 
component will tend to favor more strongly spin-polarized states. For the case of the double 
layer u = 1 systems studied by Murphy et al.,'^^ the ground state is known to already be 
an isotropic ferromagnetic state of the true spins and the addition of a parallel field would 
not, at first glance, be expected to infiuence the low energy states since they are already 
fully spin-polarized. (At a fixed Landau level filling factor, Bj_ is fixed and so the total B 
and the corresponding Zeeman energy increase with tilt). Nevertheless experiments'^^ have 
shown that these systems are very sensitive to B\\. The activation energy drops rapidly (by 
factors varying from two up to an order-of-magnitude in different samples) with increasing 
B\\ . At -Bji = -By there appears to be a phase transition to a new state whose activation gap 
is approximately independent of further increases in B\\. 

The effect of B\\ on the pseudospin system can be visualized in two different pictures, one 
microscopic, the other macroscopic. We will present the latter here. The technical details of 
the former are described elsewhere. ^^'^"^ Recent work presents a discussion of higher Landau 
levels. 

We use a gauge in which By = V x Ay where Ay = B^\ (0, 0, x). In this gauge the vector 
potential points in the z direction (perpendicular to the layers) and varies with position 
X as one moves parallel to the layers. In this gauge, the only change in the Hamiltonian 
caused by the parallel field is in the term which describes tunneling between layers. As 
an electron tunnels from one layer to the other it moves along the direction in which the 
vector potential points and so the tunneling matrix element acquires a position-dependent 
phase t ^ t e**^^ where Q = 2tt/L\\ and L\\ = (^q/ B\\d is the length associated with one flux 
quantum $o between the layers (defined in Fig.0). This modifies the tunneling Hamiltonian 
to Ht = — J d'^r h(r) ■ S(r) where h(r) 'tumbles': i.e., h(r) = 2t (cosQx, sinQx, 0). The 
effective XY model now becomes 

H = |dV|ip,|V^r-^ cos[v9(r)-gx]|, (125) 
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which is precisely the Pokrovsky-Talapov (P-T) modeP^^ and has a very rich phase dia- 
gram. For small Q and/or small Ps the phase obeys (at low temperatures) Lp{r) = Qx; the 
order parameter rotates commensurately with the pseudospin Zeeman field. However, as 
S|| is increased, the local field tumbles too rapidly and a continuous phase transition to an 
incommensurate state with broken translation symmetry occurs. This is because at large 
it costs too much exchange energy to remain commensurate and the system rapidly 
gives up the tunneling energy in order to return to a uniform state Vy? ~ which becomes 
independent of B\\ . As explained in further detail below we find that the phase transition 
occurs at zero temperature for^^'^^ 

Bl = (2£/7rci)(2i/7rp,)^/2_ ^126) 

Using the parameters for the sample of Murphy et al.^^ with weakest tunneling^^^ (Asas = 
0.45 K) and neglecting quantum fluctuation renormalizations of both t and Ps (i.e., using 
Hartree Fock) we find that the critical field for the transition is ~ 1.3T which is slightly 
larger than the observed value"^^'^^^ of 0.8 T, but still correctly corresponds to a very large 
length A graphical comparison showing the qualitative agreement between the predicted 
and observed values of the critical tilt angle for several samples is shown in Fig. (15) of the 
Chapter by Eisenstein in this volume. 

In addition to the Hartree-Fock calculations described here, we have made numerical 
exact diagonalization studies on small systems to find the critical value of the parallel field. 
Although it is difficult to extrapolate the results to the thermodynamic limit, they do confirm 
that at finite layer separation, quantum fluctuations can reduce the predicted critical parallel 
fleld.^^ 

As previously mentioned, the observed value B*^ = 0.8T corresponds in these samples to 
a large value for Ly: L\\/i ~ 20 indicating that the transition is highly collective in nature. 
We emphasize again that these very large length scales are possible in a magnetic field only 
because of the interlayer phase coherence in the system associated with condensation of a 
neutral object. 

Having argued for the existence of the commensurate-incommensurate transition, we 
must now connect it to the experimentally observed transport properties. In the commen- 
surate phase, the order parameter tumbles more and more rapidly as B\\ increases. As we 
shall see below, it is this tumbling which causes the charge gap to drop rapidly. In the 
incommensurate phase, the state of the system is approximately independent of B\\ and this 
causes the charge excitation gap to saturate at a fixed value. 

Recall that in the presence of tunneling, the cheapest charged excitation was found to be a 
pair of vortices of opposite vorticity and like charge (each having charge ±1/2) connected by 
a domain line with a constant string tension. In the absence of B\\ the energy is independent 
of the orientation of the string. The effect of is most easily studied by changing variables 
to 

^(r) = <^(r) - Qx. (127) 

This variable is a constant in the commensurate phase but not in the incommensurate phase. 
In terms of this new variable, the P-T model energy is 
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H 



cos (9) 



(128) 



We see that B\\ defines a preferred direction in tlie problem. Domain walls will want to line 
up in the y direction and contain a phase slip of a preferred sign (— 27r for Q > 0) in terms 
of the field 6. Since the extra term induced by Q represents a total derivative, the optimal 
form of the soliton solution is unchanged. However the energy per unit length of the soliton, 
which is the domain line string tension, decreases linearly with Q and hence Bf 



T = Tn 



1 - 



B\\ 



(129) 



where Tq is the tension in the absence of parallel B field given by Eq. ([l2lD , and Eff is the 
critical parallel field at which the string tension goes to zero.^^^ We thus see that by tuning 
one can conveniently control the 'chemical potential' of the domain lines. The domain 
lines condense and the phase transition occurs (in mean field theory) when the string tension 
becomes negative. 

Recall that the charge excitation gap is given by the energy of a vortex pair separated 
by the optimal distance Rq - 



/eV(4eT). From Eq. ([l22D we have that the energy gap far 
on the commensurate side of the phase transition is given by 



2E' + 



1 - 




(130) 



As B\\ increases the reduced string tension allows the Coulomb repulsion of the two vortices 
to stretch the string and lower the energy. Far on the incommensurate side of the phase 
transition, the possibility of inter layer tunneling becomes irrelevant. From the discussion of 
the previous section one can argue that the ratio of the charge gap at B\\ = to the charge 
gap at B\\ oo should be given (very roughly) by 
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a/4 



3/4 



;i3i) 



Putting in typical values of t and ps gives gap ratios in the range ~ 1.5 — 7 in qualitative 
agreement with experiment. According to the discussion of the previous section, gap ratios as 
large as ~ (tmax/i^cr)^^^ ~ 0.07(e^/e£)/ps, can be expected in the regime where the pseudospin 
texture picture applies. Here tmax is the hopping parameter at which the crossover to single- 
particle excitations occurs. Thus gap ratios as large as an order of magnitude are easily 
possible. Of course, all the discussion here neglects orbital effects (electric subband mixing) 
within each of the electron gas layers, and these will always become important at sufficiently 
strong parallel fields. [Note also that for d near d* , the system will have enhanced sensitivity 
to renormalization of parameters by electric subband mixing in tilted fields.] 
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It should be emphasized that only this highly collective picture involving large length 
scale distortions of topological defects can possibly explain the extreme sensitivity of the 
charge gap to small tilts of the B field. Recall that at Bp the tumbling length Ly is much 
larger than the particle spacing and the magnetic length. Simple estimates of the cost to 
make a local one-body type excitation (a spin-flip pair for example) shows that the energy 
decrease due to S|| is extremely small since £/L\\ is so small. Numerical exact-diagonalization 
calculations on small systems confirm the existence of this phase transition and show that 
the fermionic excitation gap drops to a much smaller value in the incommensurate phase. 

The collective excitation modes of the system in the commensurate and incommensurate 
phases appear to be quite interesting. As mentioned previously, the tumbling of the order 
parameter may produce a self-grating effect which will allow one to tune the wavevector 
which couples to optical probes. 

All of our discussion of the phase transition in a parallel field has been based on mean- 
field theory. Close to the phase transition, thermal fluctuations will be important. At flnite 
temperatures there is no strict phase transition at B*^ in the the P-T model. However there 
is a flnite temperature KT phase transition at a nearby B\i > Bp At flnite temperatures 
translation symmetry is restored^^^ in the incommensurate phase by means of dislocations 
in the domain string structure. Thus there are two separate KT transitions in this system, 
one for t = 0, the other for t 7^ and i^y > Bp Recently Read^^^ has studied this model 
at flnite temperatures in some detail and has shown that just at the critical value of B\\ 
there should be a square-root singularity in the charge gap. The existing data does not 
have the resolution to show this however. M. P. A. Fisher has pointed out^^'' that at zero- 
temperature the commensurate-incommensurate phase transition must be treated quantum 
mechanically. It is necessary to take account of the world sheets traced out by the time 
evolution of the strings which fluctuate into existence due to quantum zero-point motion. 
He has also pointed out that the inevitable random variations in the tunneling amplitude 
with position, which we have not considered at all here, cause a relevant perturbation. 



XI. SUMMARY 

We have discussed the origin of spontaneous ordering in multicomponent fractional quan- 
tum Hall effect systems. For 'real' spin at flUing factors u — 1 this spontaneous ferromag- 
netism induces a very large charge excitation gap even in the absence of a Zeeman gap. The 
charged excitations are interesting topological 'skyrmion'-like objects. Magnetic resonance 
experiments^^' have confirmed this remarkable picture which was developed analytically 
by Sondhi et al.^"^ to explain the numerical results of Rezayi.^^ 

We have discussed in detail a pseudospin analogy which shows how spontaneous inter- 
layer coherence in double layer quantum Hall systems arises. This coherent XY phase order 
occurring over long length scales is essential to explain the experimental observations of 
Murphy et al.^^ described in the Chapter in this volume by J. P. Eisenstein. The essential 
physics is condensation of a charge-neutral bosonic order parameter field (pseudospin mag- 
netization). This condensation controls the charge excitation gap and is very sensitive to 
interlayer tunneling and parallel magnetic field. 

We summarize a portion of this rich set of phenomena in the schematic zero-temperature 
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phase diagram of double layer systems shown in Fig.[0]. First consider the plane with 
^SAs = (zero tunneling). We have argued that the system develops spontaneous interlayer 
phase coherence despite the fact that the tunneling amplitude is zero. However if the 
layer spacing d exceeds a critical value d*, the system is unable to support a state with 
strong interlayer correlations and the spontaneous U{1) symmetry breaking is destroyed by 
quantum fluctuations.^^ At this same point we expect the fermionic gap A/i to collapse. 
Little is known about the nature of this quantum transition which can be viewed as arising 
from proliferation of quantum induced vortices (merons). This is similar to the quantum 
XY model, but in the present case, the merons are fractional-statistics anyons which will 
presumably change the universality class of the transition. 

At finite tunneling, the U{1) symmetry is destroyed and the quantum fluctuations are 
gapped and hence stabilized. This causes the critical layer spacing to increase as shown in 
Fig. PI . 

The third axis in the figure is the tilt of the magnetic field. Magnetic flux between the 
layers causes the order parameter to want to 'tumble'. For small tilts, the system is in a 
commensurate phase with the order parameter tumbling smoothly. However above a critical 
value of the parallel field, this tumbling costs too much exchange energy, and the system 
goes into an incommensurate phase which spontaneously breaks translation symmetry (in 
the absence of disorder). This phase transition has been observed by Murphy et al.^^'^^ 
through the rapid drop in the charge gap as the field is tilted. We have presented arguments 
that the charge gap is determined by the cost of creating a highly collective object: a pair 
of fractionally charged vortices connected by a string. It is the decrease of the string tension 
with tilt which causes the extreme sensitivity to small tilts. 

In addition to all this, there is (for zero tunneling) a finite temperature Kosterlitz- 
Thouless phase transition. If observed experimentally, this would represent the first finite 
temperature phase transition in a quantum Hall system. 

We have tried to keep the present discussion as qualitative as possible. The reader is 
directed to the many references for more detailed discussion of technical points. We close 
by noting that there are still many open questions in this very rich field concerning such 
things as edge states in multicomponent systems, proper treatment of quantum fluctuations 
for the highly collective excitations which appear to exist in these systems, and the nature 
of the phase transition at the critical value of the layer separation. 
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We are considering for the moment the case of zero temperature so that B\\ is unrenor- 
mahzed by thermal fluctuations. We also treat the problem classically. 
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TABLES 



TABLE I. Generalized Laughlin states for two component systems {S is the total spin quantum 
and * denotes a state which is not an eigenstate of Sj.). The nominal filling factors v-^ and are 
shown in parentheses for the ferromagnetic {m, m, m} states because these are are not unique (only 
their sum f is fixed). (After Ref.[22]). 



m 


m! 


n 




n 




S 


1 


1 





1 


1 


2 





1 


1 


1 


(1/2) 


(1/2) 


1 


N/2 


1 


3 





1 


1/3 


4/3 


N/4 


1 


5 





1 


1/5 


6/5 


N/3 


3 


3 





1/3 


1/3 


2/3 


* 


3 


3 


1 


1/4 


1/4 


1/2 


* 


3 


3 


1 


1/5 


1 Ir. 

1/5 


z/5 


U 


3 


3 


3 


(1/6) 


(1/6) 


1/3 


N/2 
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1/3 


1/5 


8/15 
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1/7 
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3/11 


1/11 


4/11 


Ar/4 


5 


5 
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5 
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5 


5 
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1/7 


1/7 


2/7 
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5 


5 
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1/8 
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* 


5 


5 
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1/9 


1/9 


2/9 





5 


5 
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(1/10) 


(1/10) 


1/5 


N/2 
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TABLE II. Fractional charges for some two-component fractional quantum Hall effect states, 
gives the contribution to the charge from the X' component in quasihole stateX. 
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-4/9 


1/9 


-4/9 
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TABLE III. Charged Spin Texture Energies at z^t = 1 for Double Layers Systems with Tun- 
neling, ps = ps/{e^ /tj and t = t/{e^/i) where ps is the pseudospin stiffness, t is the renormal- 
ized tunneling amplitude, £ is the magnetic length, Tq = 8ps/^ is the soliton string tension and 



^ = (2^^) is the domain wall width. 



Regime 


t < 4 X 10^ pI 


4 X 10^ pI <i< IO^Vps 


W-^/ps < t 


Nature of Charged 


Meron Pairs 


Finite Length 


Single Particle 


Excitations 




Domain Line Strings 


Excitation 


Excitation Size 


8nps 




e 


Excitation Energy 


^ 27rps 




t 
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FIGURES 



FIG. 1. Collective mode dispersion for a double-layer system at z^y = 1/2 and dji = 1.5. The 
energies of the inter-Landau- level modes are measured from huJc- The ground state is approximated 
by the (3, 3, 1) Halperin state. The plotting symbols refer to the following modes: triangles (n = 1 
sum mode); circle (n = 1 difference mode); square (n = sum mode); ; diamond (n = difference 
mode). 

FIG. 2. (a) Illustration of a simple spin-flip excitation which creates a widely separated particle 
hole pair, (b) A skyrmion spin configuration (shown in cross section). The spins gradually and 
smoothly rotate from up at the perimeter to down at the origin in a circularly symmetric spin 
textural defect. For the case of Coulomb interactions, this object costs only 1/2 the energy of the 
simple spin flip. 

FIG. 3. Knight shift measurements of the electron spin polarization of a 2DEG in the vicinity 
of filling factor 1^ = 1. (After Barrett et al., Ref.[26]). 

FIG. 4. Illustration of the path u; of a spin 5m on a unit sphere. When viewed from the origin 
of spin space, the path subtends a solid angle 0, and so the path contributes a Berry's phase SQ. 

FIG. 5. A smooth spin texture. An electron moving along the boundary of the region F in 
real space has its spin follow the path in spin space along the boundary of the region labeled uj in 

Fig.g. 

FIG. 6. Infinitesimal circuit in spin space associated with an infinitesimal circuit in real space. 

FIG. 7. Schematic conduction band edge profile for a double-layer two-dimensional electron 
gas system. Typical widths and separations are ~ d ~ lOOA and are comparable to the spacing 
between electrons within each inversion layer. 

FIG. 8. Phase diagram for the double layer QHE system (after Murphy et al.^^). Only samples 
whose parameters lie below the dashed line exhibit a quantized Hall plateau and excitation gap. 

FIG. 9. A process in double-layer two-dimensional electron gas systems which encloses flux 
from the parallel component of the magnetic field. The quantum amplitude for such paths is 
sensitive to the parallel component of the field. 

FIG. 10. The four flavors of merons. These are vortices which are right or left handed and have 
topological charge ±1/2. 
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FIG. 11. Mcron pair connected by a domain wall. Each meron carries a charge e/2 which tries 
to repel the other one. 

FIG. 12. Schematic zcro-tcmpcraturc phase diagram (with d/i increasing downwards). The 
lower surface is d* below which d > d* and the interlayer correlations are too weak to support a 
fermionic gap, Afx. The upper surface gives B*^, the commensurate-incommensurate phase bound- 
ary. As d approaches d*, quantum fluctuations soften the spin stiffness and therefore increase 
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